Nuclear spin wave quantum register for solid state quantum network nodes

ABSTRACT

A system for coupling a qubit to a register, wherein the system controls application of a protocol comprising a sequence of pulses synchronized with an RF field, the protocol further comprising a timing, a phase, and a duration of each of the pulses comprising a single qubit gate, a period and amplitude of the RF field, and a number of repeats of the sequence, so that application of the protocol controls a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment. The qubit comprises a first spin state and a second spin state both of which have a zero magnetic dipole moment; the register comprises multiple register spins having an energy level structure; and the register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing the quantum state of the qubit.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims the benefit under 35 USC 119(e) of co-pending and commonly assigned U.S. Provisional Patent Application Ser. No. 63/238,624 filed Aug. 30, 2021, by Andrei Ruskuc, Joonhee Choi, Chun-Ju Wu, and Andrei Faraon, entitled “NUCLEAR SPIN WAVE QUANTUM REGISTER FOR SOLID STATE QUANTUM NETWORK NODES,” (CIT-8694-P), which application is incorporated by reference herein.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH AND DEVELOPMENT

This invention was made with government support under Grant No(s). PHY1820790 & PHY1733907 awarded by the National Science Foundation, under Grant No(s). FA9550-18-1-0374 & FA9550-21-1-0055 awarded by the Air Force and under Grant No. N00014-19-1-2182 awarded by the Office of Naval Research. The government has certain rights in the invention.

BACKGROUND OF THE INVENTION 1. Field of the Invention

The present disclosure relates to compositions of matter useful as nuclear spin wave quantum registers and systems for implementing the same.

2. Description of Related Art

Solid-state nuclear spins surrounding individual, optically addressable qubits provide a crucial resource for quantum networks [3-6], computation [7-11] and simulation [12]. While hosts with sparse nuclear spin baths are typically chosen to mitigate qubit decoherence [13], developing coherent quantum systems in nuclear spin rich hosts enables exploration of a much broader range of materials for quantum information applications. The collective modes of these dense nuclear spin ensembles provide a natural basis for quantum storage [14]. However, utilizing them as a resource for storing quantum bits has thus far remained elusive. The present disclosure satisfies this need.

SUMMARY OF THE INVENTION

The present disclosure reports on a novel system for transferring quantum information using a qubit with zero magnetic dipole moment and indistinguishable register spins (e.g., nuclear) having an energy level structure which can be implemented in a variety of materials. The system further includes a novel protocol for controlling spin preserving interaction between the qubit and the register spins (surprisingly, despite the lack of magnetic dipole moment and the presence of noise in the system). Working embodiments described herein demonstrate the protocol can decouple the qubit from noise causing decoherence and uncontrolled/random interactions between the qubit and register, so that the spin preserving interaction can be configured to perform a variety of operations including:

-   -   Polarizing the register spins into a polarized state;     -   Generating a swap gate that transfers information between qubit         and register, and store the quantum information in the qubit in         a spin wave form described by basis states including the         polarized state and a superposition state of the register spins;         and     -   Generating a square root of swap gate used to prepare and         measure Bell states.

Devices and methods according to embodiments described herein include, but are not limited to, the following.

1. A device for coupling a qubit to a register, comprising:

a circuit for controlling application of one or more cycles of a protocol, the protocol comprising a sequence of pulses synchronized with an RF field, a timing, a phase, and a duration of each of the pulses, and a period and amplitude of the radio frequency (RF) field, wherein:

application of the protocol controls a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment;

the qubit comprises a first spin state and a second spin state both of which have a zero magnetic dipole moment;

the register comprises multiple register spins having an energy level structure,

the register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing the quantum state of the qubit, and

the pulses each comprise an electromagnetic field tuned to excite a transition between the first spin state and the second spin state.

2. The device of example 1, wherein the protocol is configured to:

suppress or cancel one or more non-exchange interactions between the register and the qubit,

suppress or cancel noise coupled to the qubit and causing decoherence of a quantum state of the qubit,

enable the coherent spin exchange interaction that performs a quantum logic gate, coherently transferring a quantum state of the qubit between the register and qubit.

3. The device of example 1 or 2, wherein the circuit controls:

application of the period of the protocol within a time period shorter than a rate of change of a magnetic noise (e.g., Overhauser field), so that the magnetic noise is quasi static during the application of the period of the protocol, the magnetic noise causing qubit decoherence and inducing a second order interaction (incoherent interaction) between the qubit and the register; and

at least one of the phase, the duration, or a time spacing of the pulses in the period so that:

-   -   one or more spin exchange interactions induced by the RF field         are preserved or maintained across the period;     -   one or more non-exchange interactions induced by the RF field         are cancelled across the period (e.g., components of the non         exchange interactions induced at different time instances in the         period cancel each other, or average to zero, over the period);     -   one or more (or any) exchange and one or more (or any)         non-exchange interactions induced by the magnetic noise are         cancelled across the period (e.g., components of these         interactions induced at different time instances in the period         cancel each other, or average to zero, over the period); and     -   the qubit decoherence induced by the magnetic noise is cancelled         over the period (e.g., decoherence induced at different time         instances in the period cancel each other, or average to zero,         over the period); and

the RF field toggling between two values of equal magnitude and opposite polarity such that:

-   -   the period is associated with a frequency of a precession of         each of the multiple register spins about a predetermined         quantization axis; and     -   the amplitude is selected for a predetermined magnitude of the         coherent spin exchange interaction between the register spins         and qubit; and

so as to form a predictable and coherent spin exchange interaction.

4. The device of any of the examples 1-3, wherein each of the single qubit gates comprises one of the pulses having the frequency and duration tuned to drive a transition between the first spin state and the second spin state.

5. The device of any of the examples claim 1-4 comprising a quantum memory, wherein the circuit:

controls application of a number of cycles the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state;

controls application of one or more of the pulses to set a quantum state of the qubit; and

controls application of a number of cycles of the protocol so as to apply a first swap gate (two qubit gate) transferring the quantum state of the qubit from the qubit to the register, thereby changing the polarized state to a corresponding state of the register spins corresponding to the quantum state; and

controls application of a number of cycles of the protocol so as to apply a second swap gate retrieving the quantum state in the qubit from the register, thereby changing the corresponding state of the register spins to the polarized state.

6. The device of any of the examples 1-5, wherein configuring the register spins in the polarized state comprises polarizing the register, which is initially in an unpolarized state comprising any configuration of excitations of the register spins, by:

(a) initializing the qubit in the first spin state by controlling application of one or more initialization pulses of an initialization electromagnetic field having a frequency tuned to initialize the quantum state of the qubit in the first spin state;

(b) applying the protocol transferring a spin excitation from the register spins to the qubit; and

(c) repeating steps (a) and (b) until all excitations of the register spins are transferred from the register to the qubit and the register spins are initialized in the polarized state, as characterized by a measurement of the qubit remaining in the first spin state after step (b).

7. The device of example 5 or 6, wherein the circuit controls application of the protocol so as to apply the first swap gate mapping (via the coherent spin exchange interaction) between the qubit and the register, such that:

if the qubit is in the first spin state, the corresponding state of the register is the polarized state,

if the qubit is in the second spin state, the corresponding state of the register is a W state, and

if the qubit is in a superposition of the first spin state and the second spin state, the corresponding state of the register is a superposition of the polarized state and the W state, and

wherein the W state is a superposition of all single spin excitation states of the register spins.

8. The device of any of the examples 1-7, wherein the circuit:

controls application of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state;

controls application of one or more of the pulses to set a quantum state of the qubit;

controls application of the protocol so as to apply a first square root of swap gate entangling the qubit with the register so as to form a Bell state; and

controls application of the protocol so as to apply a second square root of swap gate interacting with the Bell state so as to perform a measurement of the Bell state.

9. A repeater in a quantum network comprising the device of example 8.

10. A system for coupling the qubit to the register comprising the device of any of the examples 1-9, further comprising:

a photonic cavity coupled to a solid state material comprising the qubit and the register;

one or more microwave sources coupled to the qubit via a microwave waveguide, the microwave sources outputting one or more first microwave pulses and/or one or more second microwave pulses;

a radio frequency source outputting the RF field; and

one or more laser sources outputting one or more laser pulses coupled to the qubit through the photonic cavity; and wherein:

the circuit controls the one or more laser sources and the one or more microwave sources so as to:

output initialization pulses comprising at least one of the one or more laser pulses or the one or more first microwave pulses having initialization frequencies for exciting one or more transitions initializing the qubit;

apply the protocol comprising the single qubit gates comprising the second microwave pulses in synchronization with the RF field; and

output one or more readout electromagnetic fields having a readout frequency for exciting a readout transition from the second spin state to a readout state, so as to stimulate output of third pulses from the readout state.

11. The device of any of the examples 1-10, wherein:

the pulses each comprise a pi pulse or a pi/2 pulse having at least one phase selected from +x. −x., +y, or −y and

the circuit controls:

the sequence such that the period of the RF field is 2τ, and a spacing of the pulses is τ/4, and

for a given magnitude of the spin exchange interaction determined by the amplitude of the RF field, a number of repeats of the protocol that applies at least one of a swap gate transferring a quantum state between the qubit and the register, a square root of a swap gate for forming or measuring a Bell state, or that can be used to polarize the spins into a polarized state in combination with an initialization of the qubit.

12. The device of any of the examples 1-11 wherein the circuit selects the duration and the timing of each of the pulses and a toggling of the RF field to engineer the coherent spin-exchange interaction comprising:

Ŝ₊Î⁻+Ŝ⁻Î₊,

where Î₊=|↑

↓|,Î⁻=|↓

↑| are the raising and lowering operators in an effective nuclear two-level manifold of the multiple spins in the register and Ŝ₊ are similarly defined for the qubit.

13. The device of any of the examples 1-12, wherein the RF field comprises a square wave and the sequence of pulses comprise:

in a first half period τ of the square wave a sequence of the second pulses comprising:

a first pi/2 pulse having a phase +Y followed by a first pi pulse having a phase +Y, the beginning of the first pi/2 pulse and the center of the first pi pulse separated in time by τ/4;

a second pi/2 pulse immediately followed by a third pi/2 pulse, the end of the second pi/2 pulse separated in time from the center of the first pi pulse by τ/4, wherein the second pi/2 pulse has a phase −Y and the third pi/2 pulse has a phase −X;

a second pi pulse having a phase −X and following the third pi/2 pulse, a center of the second pi pulse separated in time from the center of the first pi pulse by τ/2; and

a fourth pi/2 pulse having a phase −X, wherein the end of the fourth pi/2 pulse is w separated in time from center of the second pi pulse by τ/4; and

in a second half period τ of the square wave, a repeat of the sequence of second pulses but wherein the first pi/2 pulse, the first pi pulse, and the second pi/2 pulse have opposite phase as compared to the first pi/2 pulse, the first pi pulse, and the second pi/2 pulse in the first half period, respectively.

14. The device of any of the examples 1-13, wherein the protocol de-couples the qubit from decoherence noise and random interactions caused by a nuclear Overhauser field generated by a host lattice in which the qubit is located.

15. A system for implementing a quantum register comprising the device of any of the examples 1-14 coupled to:

a spin carrying defect in a host lattice, wherein the spin carrying defect comprises the qubit and the host lattice comprises the register, or

a quantum dot in a host lattice, wherein the quantum dot comprises the qubit and the host lattice comprises the register.

16. The system of example 15, wherein the spin carrying defect is a qubit ion comprising the qubit and the register comprises a lattice of register ions surrounding the qubit ion.

17. The device of any of the examples 1-16, wherein the multiple spins in the register comprise nuclear spins and the first spin state and the second spin state comprise hyperfine electron spin states.

18. A method for coupling a qubit to a quantum register, comprising:

controlling application of a protocol comprising a sequence of pulses synchronized with an RF field, the protocol further comprising a timing, a phase, and a duration of each of the pulses comprising a single qubit gate, a period and amplitude of the RF field, and a number of repeats of the sequence, wherein:

application of the protocol controls a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment;

the qubit comprises a first spin state and a second spin state having the zero magnetic dipole moment;

the register comprises multiple register spins having an energy level structure,

the register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing the quantum state of the qubit, and

the pulses comprise an electromagnetic field tuned to excite a transition between the first spin state and the second spin state.

19. The method of claim 18, wherein the controlling further comprises:

applying the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state;

applying of one or more of the pulses to set a quantum state of the qubit;

controlling application of the protocol so as to apply a first swap gate (two qubit gate) transferring the quantum state of the qubit from the qubit to the register, thereby changing the polarized state to a corresponding state of the register spins corresponding to the quantum state; and

controlling application of the protocol so as to apply a second swap gate retrieving the quantum state in the qubit from the register, thereby changing the corresponding state of the register spins to the polarized state.

20. The method of claim 18, wherein the controlling further comprises:

controlling application of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state;

controls output of one or more of the pulses to set a quantum state of the qubit;

controls application of the protocol so as to apply a first square root of swap gate entangling the qubit with the register so as to form a Bell state; and

controls application of the protocol so as to apply a second square root of swap gate interacting with the Bell state so as to perform a measurement of the Bell state.

21. A device for controlling a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment, wherein the qubit comprises a first spin state and a second spin state having the zero magnetic dipole moment; and the register comprises multiple indistinguishable spins.

BRIEF DESCRIPTION OF THE DRAWINGS

Referring now to the drawings in which like reference numbers represent corresponding parts throughout:

FIGS. 1A-1E. Schematic of a many-body nuclear spin register for optically-coupled ¹⁷¹Yb qubits in nanophotonic cavities. FIG. 1A, Optically addressable ¹⁷¹Yb ion (yellow) surrounded by a local ensemble of nuclear spins from lattice ⁵¹V ions. The register (blue) consists of four ⁵¹V spins equidistantly spaced by 3.9 Å from the central ¹⁷¹Yb. The nuclear spin bath (grey) creates random magnetic noise termed the nuclear Overhauser field. A nanophotonic cavity enables optical initialisation and readout of the ¹⁷¹Yb ion via single-photon detection at 984 nm [15]. 675 MHz microwave pulses provide high-fidelity control of the ¹⁷¹Yb spin state. FIG. 1B, Energy level structure of ¹⁷¹Yb and ⁵¹V ions. Pulse-based control of the ¹⁷¹Yb ground-state transition (|0_(g)

↔|1_(g)

) enables engineered spin-exchange interactions with neighbouring ⁵¹V ions. The energy level structure of the spin-7/2 ⁵¹V consists of four quadratically-spaced, doubly degenerate energy levels, {|±m_(I)

}={|±1/2

, |±3/2

, |±5/2

, |±7/2

}, resulting in three distinct transitions, ω_(a,b,c)/2π=330 kHz, 660 kHz, and 991 kHz, respectively. The ω_(c) transition (dotted box) is used to implement the local nuclear spin register for quantum information storage. FIG. 1C, Effective qubit states of the nuclear spin register. The |0_(v)

and |W_(v)

states consist of all four ⁵¹V ions prepared in the |↓

=|±7/2

state and a single spin excitation equally delocalised in the |↑

=|±5/2

state, respectively. FIG. 1D, Initialisation of the nuclear spins from a thermal state into the polarised |↓↓↓↓

state. e, Transfer of a quantum state from ¹⁷¹Yb to the ⁵¹V register, storage and subsequent retrieval. Both the state initialization and transfer are enabled by robust, dynamically engineered interactions between ¹⁷¹Yb and ⁵¹V ions.

FIGS. 2A-2D. Pulse-based Hamiltonian engineering, nuclear register polarisation and spin exchange between ¹⁷¹Yb and ⁵¹V ions. FIG. 2A, Engineered spin-exchange interactions via our ZenPol sequence. Equidistant π/2 and pulses combined with a square-wave RF pulse, with magnetic field amplitude B^(RF), are periodically applied to the ¹⁷¹Yb qubit with base sequence period 2τ. FIG. 2B, ZenPol sequence spectroscopy. ¹⁷¹Yb ⁵¹V resonance is achieved for a given ⁵¹V transition (ω_(j)) when 1/τ=ω_(i)/πk with integer k. The isolated, RF-induced ω_(c) (k=5) and ω_(b) (k=3) transitions are used to polarize the multi-level nuclear spins of neighbouring ⁵¹V ions (dashed boxes). Both cases exhibit split-resonance features, attributed to the presence of two distinct ⁵¹V ensembles: the four ⁵¹V register spins (starred transitions) adjacent to the ¹⁷¹Yb qubit experience a frozen-core type detuning relative to the more distant bath. Insets: Under repeated application of the ZenPol sequence targeted at the ω_(b) or ω_(c) register transitions and interleaved with ¹⁷¹Yb initialisation, the four register spins are selectively polarised (purple lines). FIG. 2C, Spin-exchange dynamics with the four ⁵¹V register spins. The ¹⁷¹Yb qubit and ⁵¹V register spins are initialized into |1_(g)

and |0_(v)

(≡|↓↓↓↓

), respectively. Subsequently, our pulse sequence induces resonant spin exchange on the ω_(c) transition leading to oscillation between |1_(g)

|0_(v)

↔|0_(g)

|W_(v)

where |W_(v)

is a spin-wave like W-state (red markers). Inset: the rate of spin exchange scales linearly with B^(RF). With ¹⁷¹Yb in |0_(g)

, there are no spin excitations in the system and oscillations are suppressed (blue markers). A ZenPol sequence with M=10 periods and duration t_(M)=50 μs is used to realise a swap gate (black arrow). d, Spin-exchange dynamics with a single ⁵¹V nuclear spin. Three ⁵¹V spins are shelved in |±3/2

and a single spin is excited to |↑

=|±5/2

. Accordingly, the ¹⁷¹Yb qubit undergoes spin exchange with the ω_(c) transition manifold at a reduced oscillation frequency. In FIGS. 2C, 2D, solid lines are from simulations with phenomenological exponential decay constants.

FIGS. 3A-3C. Quantum information storage in the entangled nuclear spin register. FIG. 3A, Ramsey coherence time measurement. The ¹⁷¹Yb qubit is prepared in a superposition state which is subsequently swapped onto the ⁵¹V register. After waiting for a period of time, t, the superposition state is swapped back to the ¹⁷¹Yb qubit and measured in the x basis. Fast oscillations are observed at the ⁵¹Vω_(c)/2π=991 kHz frequency (inset) and the coherence is derived from the oscillation contrast. The resulting 1/e coherence decay time is measured to be 58±4 μs. Note that the wait time excludes the swap gate duration. FIG. 3B, Coherence time extension via motional narrowing of the ¹⁷¹Yb Knight field. By applying x-axis π pulses spaced by 2t_(w)=6 μs to the ¹⁷¹Yb qubit, the coherence time of the ⁵¹V register is extended to 225±9 μs.c. FIG. 3C. Further coherence enhancement via dynamical decoupling of the ⁵¹V register. In addition to the π pulses acting on ¹⁷¹Yb, two π pulses are applied to the ⁵¹V register with a variable inter-pulse delay time, 2t_(D). This rephases contributions to the detuning from the nuclear Overhauser field and leads to an extended memory time of 760±14 μs. Note that even numbers of ⁵¹Vπ pulses are necessary to return the register to the {|0 _(v)

,|W_(v)

} manifold prior to state retrieval. In FIGS. 3A-3C, solid lines are fits to Gaussian decay.

FIGS. 4A-4C. Characterization of maximally entangled ¹⁷¹Yb⁵¹V register Bell state. FIG. 4A, Parity oscillations between |Ψ⁺

and |Ψ⁻

(where |Ψ^(±)

=1/√{square root over (2)}(|1_(g)

|0_(v)

∓i|0_(g)

|W_(v)

)) revealing the Bell state coherence time. To prepare the |Ψ⁺

Bell state, a √{square root over (swap)} gate is applied to |1_(g)

|0_(v)

; subsequently during a wait time of duration t coherent parity oscillations occur between |Ψ⁺

and |Ψ⁻

at the ⁵¹Vω_(c) transition frequency. A second √{square root over (swap)} gate maps the resulting parity to ¹⁷¹Yb population. The oscillation contrast (and hence Bell state coherence) decays with a 1/e timescale of T_(2, Bell)*=8.5±0.5 μs, consistent with the ¹⁷¹Yb₂* time. FIG. 4B, Bell state coherence extension. During the parity oscillation, an XY-8 [38] decoupling sequence is applied to the ¹⁷¹Yb qubit. This leads to a significantly extended Bell state coherence time of T_(2, Bell)*=239±6 μs, limited by the ⁵¹V T₂* time measured in FIG. 3 b . FIG. 4C, Reconstructed Bell state density matrix. Diagonal entries representing populations are extracted through a sequential tomography protocol [39] (Methods). Off-diagonal matrix elements representing coherences are obtained from the parity oscillation contrast. Note that all density matrix values have been corrected to account for readout error, yielding a fidelity of 0.76±0.01. See Methods for details of the correction procedure.

FIGS. 5A-5C. Experimental setup and sequence detail. FIG. 5A, Energy level structure of ¹⁷¹Yb³⁺:YVO₄ ²F_(7/2)(0) and ²F_(5/2)(0). Initialisation into |0_(g)

involves repeated pulses on the F transition combined with consecutive pairs of π pulses applied to the A and f_(e) transitions leading to excitation into |1_(e)

. Subsequently, decay via E leads to initialisation into |0_(g)

. Optical readout relies on repeated optical π pulses on the A transition, each followed by a photon detection window during which cavity-enhanced emission via A is measured. FIG. 5B, Experimental setup. Optical control of the A and F transitions is realised via two frequency-stabilised lasers, each modulated using acousto-optic modulator (AOM) shutters. Microwave control is divided into two paths: a low frequency path consisting of 675 MHz ground state control (f_(g) transition) and RF, both generated using a single arbitrary waveform generator (AWG) channel and a high frequency path consisting of 3.4 GHz excited state microwave control (e_(e) transition). Each path is independently amplified and combined using a diplexer. The device chip and a superconducting nanowire single photon detector (SNSPD) are cooled to ≈500 mK in a cryostat. FIG. 5C, Detailed pulse sequence used for quantum state storage and retrieval. First, the ⁵¹V register and ¹⁷¹Yb qubit are initialised into |0_(v)

and |0_(g)

, respectively, as described in the main text. Subsequently, the ¹⁷¹Yb is prepared in a superposition state, via a π/2 pulse, which is swapped onto the ⁵¹V register using a ZenPol sequence resonant with the 991 kHzω_(c) ⁵¹V transition. After a wait time, t, the state is swapped back to ¹⁷¹Yb and measured in the x basis via a π/2 pulse followed by optical readout.

FIGS. 6A-6B. Randomised benchmarking and dynamical decoupling. FIG. 6A, The average fidelity of single qubit gates applied to the ¹⁷¹Yb|0_(g)

↔|1_(g)

transition is applied by application of a series of M_(gate) randomly sampled Clifford gates followed by the inverse operation (top inset). When averaged over a sufficiently large number of samples (in our case 100) it is possible to extract an average gate fidelity from the 1/e exponential decay constant, leading to f=0.9975±0.00004. FIG. 6B, We also measure the coherence time of the qubit transition using an XY-8 dynamical decoupling pulse sequence (top inset) with a fixed inter-π pulse separation of 5.6 μs and variable number of repetitions, M′. This leads to an exponential decay with 1/e time constant T₂=16±2 ms

FIGS. 7A-7C. Hartmann Hahn spectroscopy. FIG. 7A, Hartmann Hahn (HH) sequence used to perform spectroscopy of the nuclear spin environment. During the HH pulse (red), the ¹⁷¹Yb|0_(g)

↔|1_(g)

qubit transition is driven resonantly for duration t with y-phase leading to a pair of dressed states,

$\left. \left. {\left. {\left. {❘ \pm} \right\rangle = {\frac{1}{\sqrt{2}}\left( {❘0_{g}} \right.}} \right\rangle \pm {❘1_{g}}} \right\rangle \right),$

separated by energy splitting equal to the Rabi frequency, Ω. An initial −x-phase π/2 pulse prepares the ¹⁷¹Yb qubit in the |−

dressed state. When the Rabi frequency of the HH pulse is tuned to equal one of the ⁵¹V transition frequencies, the ¹⁷¹Yb is transferred into the |+

dressed state as a result of resonant population exchange (green arrows). The |+

state population is mapped to |1_(g)

with a final x-phase π/2 pulse for readout. FIG. 7B, HH spectroscopy experimental results. To identify nuclear spin resonances, both the HH pulse amplitude and duration are varied. The three evenly-spaced horizontal resonance features occurring at pulse amplitudes of 0.15, 0.3, and 0.45 (in arbitrary units, a.u.) correspond to interaction with the ω_(a), ω_(b) and ω_(c) transitions, respectively. In the no driving (Ω=0) case, the sequence probes the decoherence dynamics of the prepared |−

state i.e. it measures the Ramsey coherence time. FIG. 7C, HH spectroscopy simulation results. Simulation results agree well with the experiment, verifying that the ¹⁷¹Yb-⁵¹V interactions are dominant in our system.

FIGS. 8A-8B ZenPol sequence detail. FIG. 8A, ZenPol sequence with the toggling-frame transformation of the spin1/2{tilde over (Ŝ)}_(z) operator for the ¹⁷¹Yb qubit. The ZenPol sequence consists of a series of π and π/2 pulses about the x- and y-axes combined with a synchronously applied, square-wave RF signal with period 2τ. The Overhauser- and RF-induced interactions are determined by the toggling-frame transformations of {tilde over (Ŝ)}_(π) which are given by {tilde over (Ŝ)}_(x)f_(x) ^((OH))+{tilde over (Ŝ)}_(y)f_(y) ^((OH)) and {tilde over (Ŝ)}_(x)f_(x) ^((RF))+{tilde over (Ŝ)}_(y)f_(y) ^((RF)), respectively (see yellow and purple lines for f_(x,y) ^((OH)) and f_(x,y) ^((RF)), respectively). At the resonance condition ½τ=ω_(j)/2πk for odd integer k with ⁵¹V spin precession frequency ω_(i), the sequence realises noise-robust spin-exchange interaction with a timeaveraged Hamiltonian that only depends on the RF magnetic field amplitude. FIG. 8B, ZenPol sequence filter functions corresponding to the Fourier transforms of f_(x) ^(OH) (yellow) and f_(x) ^(RF) (purple). For a sequence with fixed τ, the peak positions determine the resonant frequencies at which ¹⁷¹Yb-⁵¹V interactions can occur. Note that the incoherent Overhauser-induced interactions occur at even-k resonances and are spectrally separated from the coherent RF-induced interactions occurring at odd-k resonances.

FIGS. 9A-9C. Polarisation of multi-level nuclear register spins. FIG. 9A, Polarisation readout by polarisation inversion (PROPI) experiments for the ⁵¹V register ω_(c) transition. The PROPI sequence performs a repeated swap operation based on the ZenPol sequence, periodically interleaved with ¹⁷¹Yb qubit readout and re-initialisation into |1_(g)

. A total of 20 polarising cycles are applied to the ω_(c) transition to polarise the ⁵¹V register into |±5/2

. As a result of register polarisation, the ¹⁷¹Yb population in |1_(g)

increases over time, indicating the accumulation of the ⁵¹V population |±5/2

(left panel). We observe that the register polarisation saturates after approximately 10 cycles. Subsequently, we perform repolarisation cycles where ¹⁷¹Yb is initialised into |0_(g)

and ⁵¹V register spins are transferred to |±7/2

with similar saturation timescale (right panel). FIG. 9B, PROPI experiments for the ⁵¹V register ω_(b) transition. Applying a ZenPol sequence resonant with the ω_(b) transition, interleaved with ¹⁷¹Yb initialisation into |1_(g)

(|0_(g)

), results in ⁵¹V register polarisation into |±5/2

(|±3/2

) as indicated by an increase (decrease) in ¹⁷¹Yb|1_(g)

population. FIG. 9C, Experimental results of ZenPol spin-exchange dynamics with varying degree of ⁵¹V register polarisation. As the number of polarisation cycles used to prepare the |0_(v)

=|±7/2

^(⊗4) state increases, the subsequent spin-exchange oscillations become more pronounced. Note that these polarisation cycles are interleaved between the ω_(b) and ω_(c) transitions.

FIGS. 10A-10D. Tunable spin-exchange rate. FIG. 10A, ZenPol sequence schematic. The square-wave RF magnetic field amplitude B^(RF) determines the ¹⁷¹Yb⁵¹V interaction strength, the pulse spacing τ/4 varies the sequence detuning from a specific ⁵¹V nuclear spin transition, and the number of ZenPol periods, M, determines the total interaction time. FIG. 10B, Simulated spin-exchange dynamics near the ω_(c) transition at k=5, probed as a function of sequence resonance frequency ω and the number of ZenPol periods, M. FIG. 10C, Measured spin-exchange dynamics showing good agreement with the numerical simulation in FIG. 10B. FIG. 10D, Experimental demonstration of tunable spin-exchange rate by varying the square-wave RF amplitude, B^(RF). When increasing B^(RF) from 0.8G to 2.0G, we observe a corresponding linear increase in the spin-exchange rate. In all cases, numerical simulations (solid lines) taking into account incomplete register polarisation, control pulse imperfections and an exponential phenomenological decay show reasonable agreement with the experimental data (markers). A simulation result without this phenomenological decay (dashed line) displays a discrepancy, which needs further investigation. See Supplementary Information for simulation details.

FIGS. 11A-11C. Direct ⁵¹V nuclear spin driving. FIG. 11A, Details of ⁵¹V nuclear spin driving scheme. To directly drive the ⁵¹V nuclear spin ω_(c) transition, a sinusoidal z-directed RP magnetic field, B_(z) ^(ocs) sin(ω_(c)t), is applied to the system at a frequency of ω_(c)/2π=991 kHz after initialising the ¹⁷¹Yb and ⁵¹V register into ═0_(g)

and |0_(v)

=|↓↓↓↓

, respectively (Drive Protocol 1). This induces an oscillating magnetic dipole moment on the ¹⁷¹Yb qubit which in turn generates an amplified transverse driving field at each ⁵¹V (Methods). Consequently, the four ⁵¹V register spins undergo independent Rabi oscillation between the |↑

=|±5/2

and |↓

=|±7/2

states. To probe the nuclear spin Rabi oscillation, the |↓

population is measured by preparing the ¹⁷¹Yb in |1_(g)

via an x-phase π pulse, performing a single swap gate and reading out the ¹⁷¹Yb population. FIG. 11B, Decoupling of magnetic field noise originating from the ¹⁷¹Yb Knight field. To improve the nuclear spin control fidelity, a train of equidistant π pulses are applied to the ¹⁷¹Yb during the driving period, thereby cancelling dephasing due to the ¹⁷¹Yb Knight field (Drive Protocol 2). Each π pulse is accompanied by a π phase shift of the sinusoidal field to ensure phase continuity of the nuclear Rabi driving and an even number of π pulses ensures the ¹⁷¹Yb state is returned to |0_(g)

at the end of the sequence (Methods). FIG. 11C, Measured ⁵¹V register Rabi oscillations using the aforementioned schemes. We observe coherent nuclear Rabi oscillations between the |↓

and |↑

states at a Rabi frequency of 2π×(7.65±0.05) kHz. An exponential decay is observed with a 1/e time constant of 280±30μ without decoupling (blue). The additional π pulses applied to the ¹⁷¹Yb qubit lead to an enhancement in control fidelity, giving a 1/e Gaussian decay time of 1040±70 μs (red). The black arrow at t≈69 μs indicates the ⁵¹V pulse used in FIG. 3 c.

FIGS. 12A-12B. ⁵¹V spin register population relaxation. FIG. 12A, Measured relaxation timescales, T₁ ^((W)), of the entangled register state, |W_(v)

, under various conditions. Top: the ⁵¹V register is prepared in the |W_(v)

state by swapping a single spin excitation from the ¹⁷¹Yb initialised into |1_(g)

. After a variable wait time, t, the ⁵¹V state is swapped back onto ¹⁷¹Yb and measured (top inset). The resulting Gaussian decay shows a 1/e relaxation time of T₁ ^((W))=39.5±1.3 μs (blue trace), limited by dephasing of the entangled |W_(v)

state. Middle: the T₁ ^((W)) lifetime can be extended by applying a series of equidistant π pulses to the ¹⁷¹Yb separated by 2t_(w)=6 μs (middle inset). This decouples the |W_(v)

state from dephasing induced by the ¹⁷¹Yb Knight field, equivalent to the coherence time extension in FIG. 12B, leading to an extended 1/e lifetime of T₁ ^((W))=127±′S (red trace). Bottom: further extension of the T₁ ^((W)) lifetime is achieved by dynamical decoupling whereby additionally two ⁵¹V π pulses are applied during the wait time with a variable pulse separation 2t_(D) (bottom inset). This gives rise to a significantly prolonged lifetime of T₁ ^((W))=640±20 μs (yellow trace), equivalent to the coherence time extension in FIG. 3 c . b, Measured relaxation timescale, T₁ ⁽⁰⁾, of the polarized register state |0_(v)

. The register is initialised in |0_(v)

and after a variable wait time, t, the ⁵¹V state is swapped onto ¹⁷¹Yb and measured (inset). We observe an exponential decay with a 1/e relaxation time of T₁ ⁽⁰⁾=054±0.08 s, likely limited by spin exchange with the bath. See Supplementary Information for detailed discussion of T₁ relaxation mechanisms.

FIGS. 13A-13D. Population measurement histograms for register fidelity characterization. FIG. 13A, Sequential tomography protocol for characterising ¹⁷¹Yb-⁵¹V populations in the basis spanned by {|0_(g)0_(v)

, |0_(g)W_(v)

, |1_(g)0_(v)

, |1_(g)W_(v)

}. Reconstructing the population probability distribution utilises Readout sequences 1 and 2, each including three consecutive ¹⁷¹Yb state readouts interleaved with single-qubit gate operations and a swap gate. FIG. 13B, Table summarizing the post-processing criteria for state attribution. Readout sequences 1 and 2 measure the {|0_(g)0_(v)

, |0_(g)W_(v)

}, and {|1_(g)0_(v)

, |1_(g)W_(v)

} populations, respectively, conditioned on the three measurement outcomes. See Methods for full details of the post-processing procedure. FIG. 13C, Reconstructed population distributions for estimating state preparation fidelity. The four basis states, {|0_(g)0_(v)

, |0_(g)W_(v)

, |1_(g)0_(v)

, |1_(g)W_(v)

}, are independently prepared by applying a combination of ¹⁷¹Ybπ pulses and swap gates to the initial |0_(g)0_(v)

state (see the insets of each subplot). Subsequently, the sequential tomography protocol (RO) is applied iteratively, alternating between Readout 1 and 2 sequences to fully reconstruct the population probability distributions. FIG. 13D, Reconstructed population distribution for the ¹⁷¹Yb⁵¹V Bell state (reproduced from FIG. 4 c ). The maximally entangled Bell state |Ψ⁺

=1/√{square root over (2)}(|1_(g)0_(v)

−i|0_(g)W_(v)

) is prepared by applying a √{square root over (swap)} gate to |1_(g)0_(v)

and measured using RO (inset). In c,d, the uncorrected and readout-corrected measurement results are presented as dashed and solid filled histograms, respectively. Populations are corrected by accounting for the swap gate error during the readout sequences (Methods).

FIGS. 14A-14C. Experimental demonstration of deterministic nuclear spin register. To demonstrate the deterministic nature of the nuclear spin register, we perform the same measurements on two additional Yb ion qubits present in the device: Ion 2 (red) and Ion 3 (yellow). Results for Ion 1 (blue) are reproduced from the main text figures for ease of comparison. a, ZenPol spectra near the ω_(c) (k=5) resonance of the ⁵¹V register spins. Notice that for all three ions, the bath and register transitions are identified at the same resonance frequencies of ω_(c) ^(bath)/2π=1028 kHz and ω_(c)/2π=991 kHz, respectively. FIG. 14B, Dynamically engineered spin-exchange dynamics between the ¹⁷¹Yb qubit and ⁵¹V register. Using constant ZenPol square-wave RF amplitude, B^(RF), we obtain equal spin-exchange rates for all three ions. FIG. 14C, Characterisation of ⁵¹V register coherence times with decoupling from the ¹⁷¹Yb Knight field. The 1/e coherence times are measured to be 225±9 μs, 273±12 μs and 261±9 μs for Ions 1, 2 and 3, respectively. All of these results demonstrate that our platform provides a nearly identical nuclear spin register for every ¹⁷¹Yb qubit in the system.

FIG. 15 . Hardware environment for implementing one or more embodiments of the present invention.

FIG. 16 . Network environment for implementing one or more embodiments of the present invention.

FIG. 17A. Flowchart illustrating a method for making a system according to one or more embodiments described herein.

FIG. 17B. Schematic showing how the protocol using Hamiltonian engineering can be used for cancelation of non-exchange contributions of the interaction

${{\hat{S}}_{z}(t)}\left( {{B_{RF}(t)} + B_{OH}} \right){\sum\limits_{register}\left\lbrack {{\alpha{{\hat{I}}_{x}^{(i)}(t)}} + {\beta{\hat{I}}_{\mathcal{z}}^{(i)}}} \right.}$

wherein interactions in green sections cancel interactions in red sections across a row.

FIG. 18 . Flowchart illustrating a method of performing operations using a protocol according to one or more embodiments.

FIG. 19 . Flowchart illustrating a method of transferring quantum information using a protocol according to one or more embodiments.

FIG. 20 . Flowchart illustrating a method of forming and measuring Bell states using a protocol according to one or more embodiments.

DETAILED DESCRIPTION OF THE INVENTION

In the following description of the preferred embodiment, reference is made to the accompanying drawings which form a part hereof, and in which is shown by way of illustration a specific embodiment in which the invention may be practiced. It is to be understood that other embodiments may be utilized, and structural changes may be made without departing from the scope of the present invention.

Technical Description

The present disclosure describes a system and method for implementing a protocol for coupling a qubit to a register. The protocol comprises a sequence of pulses synchronized with an RF field, the protocol further comprising a timing, a phase, and a duration of each of the pulses comprising a single qubit gate, a period and amplitude of the RF field, and a number of repeats of the sequence, wherein application of the protocol controls a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment. The qubit comprises a first spin state and a second spin state both of which have a zero magnetic dipole moment and the register comprises multiple register spins having an energy level structure. The register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing the quantum state of the qubit. The system further typically includes a source of the pulses comprising an electromagnetic field tuned to excite a transition between the first spin state and the second spin state.

The quantum memory can be implemented using nuclear spin-wave like states that can be implemented in a variety of (e.g., solid state) material systems. In typical examples, the quantum register is implemented using utilizing high spin, spectrally-indistinguishable, dense, lattice nuclear spins surrounding solid-state qubits. The control protocols induce coherent interaction between a central solid-state qubit and surrounding lattice ion nuclear spins. Specifically, the protocols are used to generate entangled states between the solid-state qubit and local nuclear ensemble and to implement a deterministic quantum register using the same ensemble. These features are vital ingredients for building large-scale multi-node quantum networks.

The following examples demonstrate an embodiment of the auxiliary nuclear-spin-based quantum register using single rare-earth ion qubits, although other material systems (including non-nuclear spin systems) may be used.

1. First Example: System Implemented in Yb:YVO

FIGS. 1A-1E illustrate a highly coherent, optically addressed ¹⁷¹Yb³⁺ qubit doped into a nuclear spin-rich yttrium orthovanadate crystal [15] combined with a robust quantum control protocol to manipulate the multi-level nuclear spin states of neighbouring ⁵¹V⁵⁺ lattice ions. Via a dynamically-engineered spin exchange interaction, this nuclear spin ensemble is polarized to generate collective spin excitations, and subsequently used to implement a long-lived quantum memory. Unlike conventional, disordered nuclear spin based quantum memories [16-24], the platform is deterministic and reproducible, ensuring identical quantum registers for all ¹⁷¹Yb³⁺ qubits. The approach provides a framework for utilising the complex structure of dense nuclear spin baths, paving the way for building large-scale quantum networks using single rare-earth ion qubits [15,25-28].

The hyperfine levels of single ¹⁷¹Yb³⁺ ions doped into yttrium orthovanadate (YVO₄), coupled to nanophotonic cavities, form high-quality optically addressable qubits [15]. The surrounding ⁵¹V⁵⁺ lattice ion nuclear spins generate a noisy magnetic field environment due to their large magnetic moment and high spin (I=7/2). Coherent ¹⁷¹Yb qubit operation is enabled by magnetically-insensitive transitions, leading to long coherence times (16 ms) and high gate fidelities (0.99975) (FIG. 6 ). Whilst decoupling from sources of magnetic noise achieves an excellent operating regime for the ¹⁷¹Yb qubit, the ⁵¹V nuclear spins also provide a readily accessible, local resource for quantum information storage due to their inherently weak interactions with the environment. To date, most research regarding host nuclear spin utilisation has focused on several spectrally distinguishable impurity nuclear spins coupled to a localised electronic spin, e.g. ¹³C coupled to colour centres in diamond or ²⁹Si coupled to defects in silicon carbide, rare-earth ions, quantum dots or donor qubits in silicon [10,16-24]. Recently, a regime consisting of a large number of indistinguishable nuclear spins coupled to the delocalised electronic spin in a quantum dot has also been explored [29,30]. In contrast, the system described herein addresses a new regime where a small, deterministic cluster of spectrally indistinguishable nuclear spins are coupled to a single localized electronic spin. Specifically, the ¹⁷¹Yb electronic wavefunction is confined to the lattice site, and the YVO₄ crystal consists of highly isotopically pure, 99.8% ⁵¹V, nuclear spins. This confined, dense nuclear spin ensemble could be used as a deterministic local quantum processor by creating and manipulating entangled states, such as collective spin wave-like excitations, for near-term quantum applications. Critically, interfacing with these nuclear spins whilst preserving high qubit coherence necessitates the development of novel quantum control protocols using magnetically insensitive transitions that are robust against environmental noise.

At zero-magnetic field the ¹⁷¹Yb ground state contains a pair of levels |0_(g)

and |1_(g)

, separated by 675 MHz, which form our qubit [31] (FIG. 1 b ). The |1_(g)

population is optically read out via a series of π pulses at 984 nm, each followed by time-resolved detection of resonant photon emission (FIG. 5 ). This is enabled by coupling the ¹⁷¹Yb ion to a nanophotonic cavity leading to high transition cyclicity, reduced optical lifetime and high photon collection efficiency [15]. The local crystalline environment consists of ⁸⁹Y, ⁵¹V and ¹⁶O ions. Of these, ⁵¹V with nuclear spin 7/2 has the largest magnetic dipole moment and zero-field structure due to a quadrupole interaction with the lattice electric field [32]. This leads to four quadratically-spaced, doubly degenerate energy levels, {|±m₁

}={|±1/2

, |±5/2

, |±5/2

, |±7/2

}, and three magnetic-dipole allowed transitions between these levels ω_(a),ω_(b),ω_(c) (FIG. 1 b ).

Local ⁵¹V ions are categorised into two complementary ensembles: the register and the bath. The register spins fulfil two conditions: (1) they are constituents of the frozen core: a set of ⁵¹V ions spectrally distinguished from the bath due to proximity to ¹⁷¹Yb; (2) the ¹⁷¹Yb⁵¹V interaction Hamiltonian can drive transitions between their quadrupole levels. As shown later, experimental evidence suggests that the register consists of four ⁵¹V spins, equidistant from the central ¹⁷¹Yb. At zero field, the ¹⁷¹Yb|0_(g)

, |1_(g)

states have no magnetic dipole moment and thus interactions with ⁵¹V register spins are forbidden to first order. However, a weak ¹⁷¹Yb dipole moment is induced by a random magnetic field originating from the bath (the nuclear Overhauser field, with z component B_(π) ^(OH)), giving rise to an effective ¹⁷¹Yb-⁵¹V register interaction. Specifically, a second-order pertur-bation analysis yields the following Hamiltonian:

${{\hat{H}}_{int} = {{\overset{\hat{\sim}}{S}}_{z}B_{z}^{OH}{\sum\limits_{i \in {register}}\left( {{a_{x}{\hat{I}}_{x}^{(i)}} + {a_{z}{\hat{I}}_{z}^{(i)}}} \right)}}},$

where {tilde over (Ŝ)}_(π) is the ¹⁷¹Yb qubit operator along the axis in a weakly perturbed basis, Î_(x,z)(i) are the nuclear spin-7/2 operators along the x,z axes, and a_(x,z) are the coupling coefficients (Supplementary Information). Note that B_(z) ^(OH) varies randomly in time as the bath changes state in a stochastic fashion, rendering this interaction Hamiltonian unreliable for register quantum state manipulation and requiring an alternative approach. To this end, we develop a protocol to generate a deterministic ¹⁷¹Yb-⁵¹V interaction via Hamiltonian engineering, which will be elaborated later.

An additional challenge is presented by the spectral indistinguishability of the register spins, necessitating storage in collective states. As originally proposed for quantum dots [14], single spin excitations of a polarised nuclear spin ensemble can be used for quantum information storage. These states are often termed spin waves or nuclear magnons and are generated by spin-preserving exchange dynamics. Preparing these collective nuclear spin states relies firstly on initialising the thermal register ensemble into a pure state, |0_(v)

=|↓↓↓↓

, where {|↑, |↓

}={|±5/2

, |±7/2

} is a two-level sub-manifold of the nuclear spin-7/2 ⁵¹V ion (FIG. 1 c, d ). Next, with access to exchange dynamics and ¹⁷¹Yb initialised |1_(a)

, we can transfer a single excitation from the ¹⁷¹Yb to the register. It is noted that the excitation is delocalised equally across the four register spins due to coupling homogeneity as determined by the lattice geometry, thus naturally realising the entangled four-body W-state |W_(v)

[33] given by

$\left. {❘W_{v}} \right\rangle = \frac{\left. {\left. {\left. {{\left. {❘\left. \uparrow\downarrow\downarrow\downarrow \right.} \right\rangle +}❘\left. \downarrow\uparrow\downarrow\downarrow \right.} \right\rangle + {❘\left. \downarrow\downarrow\uparrow\downarrow \right.}} \right\rangle + {❘\left. \downarrow\downarrow\downarrow\uparrow \right.}} \right\rangle}{2}$

(FIG. 1C). If the ¹⁷¹Yb qubit is initialised into |0_(g)

there are no spin excitations in the system and the ⁵¹V register remains in |0_(v)

. Crucially, these dynamics realise a quantum swap gate between a target state prepared by the ¹⁷¹Yb qubit, |ψ

=α|0_(g)

+β|1_(g)

, and the |0_(v)

state of the ⁵¹V register, leading to

(α|0_(g)

+β|1_(g)

)|0_(v)

→|0_(g)

(α|0_(v)

+β|W _(v)

).

After waiting for a certain period of time, the stored quantum state can be retrieved from the nuclear register by applying a second swap gate (FIG. 1E). Note that the spin-wave like state |W_(v)

of the nuclear ensemble is being utilized as a constituent of the quantum memory basis.

To realise this storage protocol, the ¹⁷¹Yb⁵¹V spin-exchange interactions are rendered independent from the random, bath-induced dipole moment (equation (1)). Established pulse-based methods used to generate such interactions, e.g. Hartmann Hahn [34] and PulsePol [35], do not suit the requirements of the present application as they are susceptible to random noise from the bath (FIG. 7 ). To this end, the present invention uses a framework for robust dynamic Hamiltonian engineering [36] to design a new sequence tailored for qubits with no intrinsic magnetic moment (subsequently referred to as ZenPol for ‘zero first order Zeeman nuclear-spin polarisation’). ZenPol comprises equidistant π/2 and π pulses combined with a synchronous, z-directed, square-wave RF magnetic field with tuneable amplitude, B^(RF), and period 2τ (FIG. 2 a ). The sequence is repeated M times leading to a total interrogation duration of t_(M)=2τM. The RF field induces an alternating ¹⁷¹Yb magnetic dipole moment, thereby generating a similar ¹⁷¹Yb-⁵¹V interaction as B_(π) ^(OH) in equation (1) but in a controlled manner. The sequence is synchronised with the ⁵¹V precession at one of the nuclear spin transition frequencies, ω_(i), by satisfying

${\frac{1}{2\tau} = \frac{\omega_{j}}{2\pi k}},$

with k an odd integer (FIG. 8 ). At this resonance condition the leading-order dynamics are understood by considering the temporal interference between time-varying ¹⁷¹Yb spin operators and ⁵¹V precession in the interaction picture (Methods). The ZenPol sequence is designed such that RF-induced spin-preserving dynamics interfere constructively, while all other dynamics, including the bath-induced incoherent interactions, undergo destructive interference. As a result, the ¹⁷¹Yb-⁵¹V interaction is governed by the following timeaveraged effective Hamiltonian

${{\hat{H}}_{avg}b_{({k,\omega_{j}})}B^{RF}{\sum\limits_{i \in {register}}\left( {{{\overset{\hat{\sim}}{S}}_{+}{\overset{\hat{\sim}}{I}}_{-}^{(i)}} + {{\overset{\hat{\sim}}{S}}_{-}{\overset{\hat{\sim}}{I}}_{+}^{(i)}}} \right)}},$

where b_((k,ω) _(j) ₎ is a k-dependent prefactor for the ω_(i) transition, Î₊=|↑

↓|, Î⁻=|↓

(↑| are the raising and lowering operators in an effective nuclear two-level manifold and {tilde over (Ŝ)}₊ are similarly defined for the ¹⁷¹Yb qubit (Methods). While the nuclear spin can stochastically occupy either the {|+m_(I)

} or {|−m_(I)

} manifold of states, the protocol described herein is insensitive to this sign. Moreover, this pulse sequence operates at zero magnetic field where a long ¹⁷¹Yb coherence time can be maintained; it is insensitive to the presence of random noise from the bath; and is also robust to experimental imperfections, e.g. pulse rotation errors.

2. Example Protocol for the First Example

The ZenPol sequence is used to perform spectroscopy of the ¹⁷¹Yb nuclear spin environment. FIG. 2 b shows a ZenPol spectrum obtained by initialising the ¹⁷¹Yb into |0_(g)

, applying an M=30 period ZenPol sequence with variable inter-pulse spacing (τ/4) and reading out the ¹⁷¹Yb population. As a result of the engineered exchange interaction, the |0_(g)

population decreases significantly at expected τ values corresponding to the odd-k ⁵¹V resonances (red line, FIG. 2 b ). Even-k resonances are also observed even in the absence of the RF field, which are attributed to the incoherent interaction dominated by the random nuclear Overhauser field (blue line, FIG. 2 b ).

In particular, all the odd-k resonances are split near each isolated ⁵¹V transition (dotted boxes, FIG. 2 b ). For example, resonance frequencies of {660 kHz, 685 kHz} and {991 kHz, 1028 kHz} are identified around the ω_(b) (k=3) and ω_(c) (k=5) transitions, respectively. In both cases, the higher-frequency resonance agrees well with literature values extracted from NMR on YVO crystals (685 kHz, 1027 kHz) [32]. The presence of two nuclear spin ensembles is postulated: a distant large ensemble with unperturbed frequency (constituents of the bath) and a local small ensemble with a frequency shift due to crystalline strain in the vicinity of the ¹⁷¹Yb ion (the register).

Polarisation of the entire nuclear spin register relies on repeated application of the ZenPol sequence, resonant with a targeted transition, interleaved with reinitialisation of the ¹⁷¹Yb qubit leading to unidirectional transfer of ⁵¹V population. (FIG. 5 c ). Since a the spin-7/2 ⁵¹V ions have four energy levels, high fidelity initialisation by independently polarising different transitions is achieved with different values of τ. For example, to prepare the register spins in |±7/2

=|↓

, a pair of ZenPol sequences is repeatedly applied which first polarise into |±5/2

using the ω_(b) transition, and then subsequently into ±|7/2

using the ω_(c) transition (FIG. 9 ). The data confirms that both ω_(b) and ω_(c) transitions of the ⁵¹V register are successfully polarised as indicated by the near-complete disappearance of the initial resonances (insets, FIG. 2 b ). Note that the resonances at 685 kHz and 1028 kHz are unaffected, corroborating the existence of two distinct ⁵¹V ensembles discussed above. The ω_(a) transition is not directly addressed by the ZenPol sequence due to spectral overlap with other resonances, however, this does not limit our polarisation fidelity, estimated to be ≈84%, as discussed in Supplementary Examples.

After initialising all four register ⁵¹V spins into a polarized state |0_(v)

=|↓↓↓↓

, the ZenPol sequence can also induce coherent oscillations of a single spin excitation between the ¹⁷¹Yb ion and the polarised ⁵¹V ensemble. FIG. 2 c shows the ¹⁷¹Yb population as a function of sequence period, M, when the single-spin exchange is targeted at the ω_(c) transition. With ¹⁷¹Yb initialised in |1_(g)

, the quantum state evolves according to:

❘ψ(t_(m))⟩ = ❘1_(g)⟩❘0_(v)⟩cos (J_(ex)t_(m)/2) − i❘0_(g)⟩❘W_(v)⟩sin (J_(ex)t_(M)/2)

with spin-exchange rate J_(ex)=4b_((5,ω) _(c) ₎B^(RF) (red, FIG. 2 c ). Note that when J_(ex)t_(M)=π, the sequence realises a swap gate (black arrow, FIG. 2 c ), whereby a single-spin excitation is completely transferred to the register, i.e., |1_(g)

|0_(v)

→|0_(g)

|W_(v)

. Furthermore, a can be accurately controlled by var in B^(RF), allowing for fidelity optimisation of the swap gate (inset, FIG. 2 c ). By contrast, with ¹⁷¹Yb initialised in |0_(g)

, exchange interactions are forbidden and thus oscillations are suppressed (blue, FIG. 2 c ).

The spin-exchange rate is collectively enhanced by a factor of √{square root over (N)}, where N is the number of indistinguishable spins forming the register. This is verified by controlling the number of spins in the ω_(c) transition manifold and measuring the effect on J_(ex). This is implemented by first emptying the ω_(c) manifold via the application of downward-polarising ZenPol sequences, thereby pumping all four spins to |±3/2

and |±1/2

. Subsequently, a single excitation is performed on the ω_(b) transition to flip one spin from |±3/2

to |↑

(=|±5/2

), leading to N=1 spins in the ω_(c) manifold. Applying a ZenPol sequence resonant with the ω_(c) transition, it is found that the resulting exchange frequency is reduced by a factor of ≈√{square root over (4)} (FIG. 2 d ); according to the YVO₄ lattice structure, the register likely consists of the second-nearest shell of four equidistant ⁵¹V ions (Supplementary Examples). This assumption is supported by close agreement between experiment and numerical simulation in all cases (FIG. 10 ).

3. Example Implementation of the First Example as Quantum Memory

To evaluate the performance of the ⁵¹V register as a quantum memory, its information storage times are characterized under various conditions. Specifically, a superposition state is first transferred from the ¹⁷¹Yb qubit,

$\left. \left. {\left. {\frac{1}{\sqrt{2}}\left( {❘0_{g}} \right.} \right\rangle + {i{❘1_{g}}}} \right\rangle \right),$

to the ⁵¹V register via the ZenPolbased swap gate. Subsequently, the transferred state

${\frac{1}{\sqrt{2}}\left( {❘0_{v}} \right)} + {\left. ❘W_{v} \right\rangle)}$

is stored for a variable wait time, t, before being swapped back to the ¹⁷¹Yb and measured along the x-axis, thereby probing the coherence of the final state. As shown in FIG. 3 a , there is a sinusoidal oscillation of the ¹⁷¹Yb population, modulated by a Gaussian coherence decay, whose contrast vanishes with a 1/e time of T₂*=58±4 μs. This oscillation has a frequency of ω_(c)/2π=991 kHz, originating from relative phase accumulation between |0_(v)

and |W_(v)

during the wait time. The coherence time of the ⁵¹V register is predominantly limited by local magnetic field noise from two sources: a fluctuating ¹⁷¹Yb dipole moment (¹⁷¹Yb Knight field) and the nuclear Overhauser field (Supplementary Information). As shown in FIG. 3 b , the noise created by ¹⁷¹Yb can be effectively decoupled from the register by periodically flipping the ¹⁷¹Yb magnetic dipole orientation via a series of π pulses. Similar to the motional narrowing effect [37], the neutralization of the dipole moment arrests undesired phase diffusion of the register, leading to an increased 1/e coherence time of T₂*=225±9 μs. The coherence time can be further extended by performing dynamical decoupling on the ⁵¹V register to mitigate the decoherence effect of the nuclear spin bath. This relies on applying ⁵¹Vπ pulses resonant with the ω_(c) transition whilst leaving the bath unperturbed (FIG. 11 ). In FIG. 3 c , we apply two ⁵¹Vπ pulses with variable inter-pulse delay, combined with periodic pulses applied to the ¹⁷¹Yb qubit, significantly extending the the 1/e coherence time to T₂=760±1.4 μs.

The population relaxation times of the |0_(v)

and |W_(v)

states are characterized with measured lifetimes of T₁ ⁽⁰⁾=0.54±0.08 s and T₁ ^((W)):=39.5±−1.3 μs, respectively. Due to the entangled nature of the |W_(v)

state, T₁ ^((W)) is limited by dephasing and is extended to 127±8 μs and 640±20 μs by applying the same decoupling sequences as in FIG. 3 b,c respectively (FIG. 12 ). These dephasing processes can be sensitive to the stochastic occupation of the |+m_(l)

and |−m_(l)

states, depending on the degree of noise correlation between the four register spins (See Supplementary Examples).

4. Example Bell State Generation Using the First Example

The multi-spin register is benchmarked by characterizing fidelities of ¹⁷¹Yb ⁵¹V Bell state generation and detection, serving as a vital component of the quantum repeater protocol [3]. In particular, the maximally entangled Bell state

$\left. ❘\Psi^{+} \right\rangle = {{\frac{1}{\sqrt{2}}\left( 1_{g} \right\rangle\left. ❘\text{?} \right\rangle} - {\left. i❘0_{g} \right\rangle\left. ❘\text{?} \right\rangle)}}$ ?indicates text missing or illegible when filed

can be prepared by initialising the system in |1_(g)

|0_(v)

and applying a √{square root over (swap)} gate based on the ZenPol sequence satisfying J_(ex)t_(M)=π/2 (equation (6)). The Bell state coherence is evaluated by monitoring the contrast of oscillation between a given Bell state and its parity conjugate [40]. In our system, the free evolution of |Ψ⁺

gives rise to a parity oscillation at frequency ω_(c) with

$\left. \left. {\left. {\left. {\left. {\left. {❘\Psi^{-}} \right\rangle = {\frac{1}{\sqrt{2}}\left( {❘1_{g}} \right.}} \right\rangle{❘0_{v}}} \right\rangle + {i{❘0_{g}}}} \right\rangle{❘W_{v}}} \right\rangle \right)$

(See Supplementary Examples). This oscillation is read out by applying a second √{square root over (swap)} gate to the system, encoding the parity into ¹⁷¹Yb population. FIG. 4 shows the measured parity oscillations decaying with a 1/ε time of T_(2,Bell)*=8.5±0.5 μs, limited by the T₂* dephasing time of the ¹⁷¹Yb qubit [15]. To improve the coherence, an XY−8 decoupling sequence [38] is applied to the ¹⁷¹Yb, leading to an enhanced value of T_(2,Bell)*=239±6 μs (FIG. 4 b ); this timescale is similar to that in FIG. 3 b , indicating that the Bell state coherence is likely limited by the T₂* dephasing time of the ⁵¹V register.

In order to estimate the Bell state preparation fidelity, defined as

=(Ψ⁺|ρ|Ψ⁺), a sequential tomography protocol [39] is performed to reconstruct the system density matrix ρ in the effective manifold spanned by four states {|0_(g)0_(v)

, |0_(g)W_(v)

, |1_(g)0_(v)

, 1_(g)W_(v)

} (FIG. 13 ). Taking into account errors in state readout, a corrected Bell state fidelity of 0.76±0.01 is obtained, as summarized in FIG. 4 c (the uncorrected fidelity is measured to be 0.61±0.01). Without being bound by a scientific theory, this may be limited by a combination of incomplete register initialisation, imperfect Hamiltonian engineering and detrimental dephasing of the register during Bell state generation. See Methods and Supplementary Examples for detailed discussions including error analysis.

Thus, the above described examples demonstrate a noise-robust control protocol to coherently manipulate the local ⁵¹V nuclear ensemble surrounding a single optically-addressed ¹⁷¹Yb spin, enabling the polarisation of the high spin

$\left( {I = \frac{\text{?}}{2}} \right)$ ?indicates text missing or illegible when filed

nuclear register, the creation of collective spinwave excitations, and the preparation of maximally entangled Bell states. Based on these capabilities, it is shown that the local nuclear spins realise an ensemble-based quantum memory exhibiting long coherence times. Crucially, this memory is deterministic and reproducible in that every ¹⁷¹Yb ion doped into a YVO₄ crystal accesses a near-identical nuclear register in its local environment (FIG. 14 ). This resource will enable the implementation of multi-node quantum network architectures using rare-earth ions with both enhanced connectivity and large-scale entanglement [3]. Furthermore, realising coherent quantum systems using dense lattice nuclear spins demonstrates the possibility of implementation in other materials for quantum information applications [13]. These multi-level nuclear spin ensembles offer an attractive, highly controllable platform to investigate the many-body dynamics of a much larger Hilbert space, paving the way for application of solidstate, noisy intermediate-scale quantum (NISQ) devices in the context of quantum simulation [12,30]

REFERENCES FOR EXAMPLES 1-4

The following references are incorporated by reference herein

-   [1] Awschalom, D. D., Hanson, R., Wrachtrup, J. & Zhou, B. B.     Quantum technologies with optically interfaced solid-state spins.     Nat. Photonics 12, 516-527 (2018). -   [2] Chatterjee, A. et al. Semiconductor qubits in practice. Nat.     Rev. Phys. 3, 157-177 (2021). -   [3] Briegel, H. J., Dür, W., Cirac, J. I. & Zoller, P. Quantum     repeaters: The role of imperfect local operations in quantum     communication. Phys. Rev. Lett. 81, 5932-5935 (1998) -   [4] Hensen, B. et al. Loophole-free Bell inequality violation using     electron spins separated by 1.3 kilometres. Nature 526, 682-686     (2015). -   [5] Bhaskar, M. K. et al. Experimental demonstration of     memory-enhanced quantum communication. Nature 580, 60-64 (2020). -   [6] Pompili, M. et al. Realization of a multi-node quantum network     of remote solid-state qubits. Science 372, 259 264 (2021). -   [7] Waldherr, G. et al. Quantum error correction in a solidstate     hybrid spin register. Nature 506, 204-207 (2014). -   [8] Taminiau, T. H., Cramer, J., Van Der Sar, T., Dobrovitski, V. V.     & Hanson, R. Universal control and error correction in multi-qubit     spin registers in diamond. Nat. Nanotechnol. 9, 171-176 (2014). -   [9] Zhong, M., Ahlefeldt, R. L. & Sellars, M. J. Quantum information     processing using frozen core Y3+ spins in Eu3+:Y2SiO5. New J. Phys.     21 (2019). -   [10] Bradley, C. E. et al. A Ten-Qubit Solid-State Spin Register     with Quantum Memory up to One Minute. Phys. Rev. X. 9, 31045 (2019). -   [11] Kinos, A. et al. Roadmap for Rare-earth Quantum Computing.     Preprint at https://arxiv.org/abs/2103.15743 (2021) -   [12] Randall, J. et al. Observation of a manybody-localized discrete     time crystal with a programmable spin-based quantum simulator.     Preprint at http://arxiv.org/abs/2107.00736 (2021). -   [13] Wolfowicz, G. et al. Quantum guidelines for solid-state spin     defects. Nature Reviews Materials (2021). -   [14] Taylor, J. M., Marcus, C. M. & Lukin, M. D. Long-Lived Memory     for Mesoscopic Quantum Bits. Phys. Rev. Lett. 90, 4 (2003). -   [15] Kindem, J. M. et al. Control and single-shot readout of an ion     embedded in a nanophotonic cavity. Nature 580, 201-204 (2020). -   [16] Gurudev Dutt, M. V. et al. Quantum Register Based on Individual     Electronic and Nuclear Spin Qubits in Diamond. Science 316,     1312-1317 (2007). -   [17] Kolkowitz, S., Unterreithmeier, Q. P., Bennett, S. D. &     Lukin, M. D. Sensing distant nuclear spins with a single electron     spin. Phys. Rev. Lett. 109, 1-5 (2012). -   [18] Taminiau, T. H. et al. Detection and control of individual     nuclear spins using a weakly coupled electron spin. Phys. Rev. Lett.     109, 137602 (2012). -   [19] Zhao, N. et al. Sensing single remote nuclear spins. Nat.     Nanotechnol. 7, 657-662 (2012). -   [20] Metsch, M. H. et al. Initialization and Readout of Nuclear     Spins via negatively charged Silicon-Vacancy Center in Diamond.     Phys. Rev. Lett. 122, 190503 (2019). -   [21] Bourassa, A. et al. Entanglement and control of single nuclear     spins in isotopically engineered silicon carbide. Nat. Mater. 19,     1319-1325 (2020). -   [22] Hensen, B. et al. A silicon quantum-dot-coupled nuclear spin     qubit. Nat. Nanotechnol. 15, 13-17 (2020). -   [23] Kornher, T. et al. Sensing Individual Nuclear Spins with a     Single Rare-Earth Electron Spin. Phys. Rev. Lett. 124, 170402(2020) -   [24] Wolfowicz, G. et al. 29Si nuclear spins as a resource for donor     spin qubits in silicon. New J. Phys. 18 (2016). -   [25] Utikal, T. et al. Spectroscopic detection and state preparation     of a single praseodymium ion in a crystal. Nat. Commun. 5, 1-8     (2014). -   [26] Siyushev, P. et al. Coherent properties of single rareearth     spin qubits. Nat. Commun. 5, 1-6 (2014). -   [27] Zhong, T. et al. Optically Addressing Single Rare-Earth Ions in     a Nanophotonic Cavity. Phys. Rev. Lett. 121, 183603 (2018) -   [28] Chen, S., Raha, M., Phenicie, C. M., Ourari, S. &     Thompson, J. D. Parallel single-shot measurement and coherent     control of solid-state spins below the diffraction limit. Science     370, 592-595 (2620). -   [29] Gangloff, D. A. et al. Quantum interface of an electron and a     nuclear ensemble. Science 364, 52-66 (2019). 1812.07540, -   [30] Gangloff, D. A. et al. Revealing beyond-mean-field correlations     in a nuclear ensemble via a proxy qubit. Preprint at     https://arxiv.org/abs/2012.11279 (2020). -   [31] Kindem, J. M. et al. Characterization of Yb3+171:YVO4 for     photonic quantum technologies. Phys. Rev. B 98, 1-10 (2018). -   [32] Bleaney, B., Gregg, J. F., De Oliveira, A. C. & Wells, M. R.     Nuclear magnetic resonance of 51 V(I=7/2): in lanthanide     vanadates: II. The nuclear electric quadrupole interaction. J.     phys., C, Solid state phys. 15, 5293-5303 (1982). -   [33] Weimer, H., Yao, N. Y. & Lukin, M. D. Collectively enhanced     interactions in solid-state spin qubits. Phys. Rev. Lett. 110, 1-5     (2013). -   [34] Hartmann, S. R. & Hahn, E. L. Nuclear double resonance in the     rotating frame. Phys. Rev. 128, 2042-2053 (1962). -   [35] Schwartz, I. et al. Robust optical polarization of nuclear spin     baths using Hamiltonian engineering of nitrogen vacancy center     quantum dynamics. Sci. Adv. 4, 1-8 (2018). -   [36] Choi, J. et al. Robust Dynamic Hamiltonian Engineering of     Many-Body Spin Systems. Phys. Rev. X. 10, 31002 (2019) -   [37] Bauch, E. et al. Ultralong Dephasing Times in Solid State Spin     Ensembles via Quantum Control. Phys. Rev. X. 8, 031025 (2018). -   [38] Gullion, T., Baker, D. B. & Conradi, M. S. New, compensated     Can-Purcell sequences. J. Magn. Reson. 89, 479-484 (1990). -   [39] Kalb, N. et al. Entanglement distillation between solidstate     quantum network nodes. Science 356, 928-932 (2017). -   [40] Levine, H. et al. High-Fidelity Control and Entanglement of     Rydberg-Atom Qubits. Phys. Rev. Lett. 121, 1-6(2018) -   [41] Zhong, T., Rochman, J., Kindem, J. M., Miyazono, E. &     Faraon, A. High quality factor nanophotonic resonators in bulk     rare-earth doped crystals. Opt. Express 24, 536, (2016). -   [42] Zhong, T. et al. Nanophotonic rare-earth quantum memory with     optically controlled retrieval. Science 1395, 1392-1395 (2017). -   [43] Dreyer, R. W. P. et al. Laser Phase and Frequency Stabilization     Using an Optical Resonator. Appl. Phys. B 31, 97-105 (1983). -   [44] Slichter, C. P. Principles of Magnetic Resonance     (Springer-Verlag, New York, 1992), 3rd edn. -   [45] Degen, C. L., Reinhard, F. & Cappellaro, P. Quantum sensing.     Rev. Mod. Phys. 89, 1-39 (2017). -   [46] Bernien, H. et al. Heralded entanglement between solidstate     qubits separated by three metres. Nature 497, 86-90 (2013). -   [47] Nguyen, C. T. et al. An integrated nanophotonic quantum     register based on silicon-vacancy spins in diamond. Phys. Rev. B     100, 1-19 (2019).

5. Supplementary Example Methods for Implementing the First Example

a. Experimental Setup

FIG. 5 illustrates a schematic of the complete experimental setup used for characterization of the first example.

The YVO₄ crystal used in this project was cut and polished from an undoped boule (Gamdan Optics) with a residual total ¹⁷¹Yb concentration of 140 ppb. Nanophotonic cavities were fabricated from this material using focused ion beam milling, see [41,42] for more detail on this process. The cavity used in this work has a Qfactor of ≈10,000 leading to Purcell enhancement and consequent reduction of the ¹⁷¹Yb excited state lifetime from 267μ to 23μ as described and measured in [15] and >99% of ion emission coupling to the cavity mode. The reduced optical lifetime enables detection of single ¹⁷¹Yb ions. The cavity is undercoupled with κ_(in)/κ≈0.14 leading to 14% of emitted light entering the waveguide mode. Waveguide-freespace coupling is achieved via angled couplers with an efficiency of ≈25% and the end-to-end system efficiency (probability of detecting an emitted photon) is ≈1%.

The device sits on the still-plate of a ³He cryostat (Bluefors LD-He250) with base temperature of 460 mK. Optical signals are fed into the fridge through optical fibre and focused onto the device with an aspheric lens doublet mounted on a stack of x-y-z piezo nanopositioners (Attocube). The device is tuned on-resonance with the ¹⁷¹Yb optical transitions via nitrogen condensation. Residual magnetic fields are cancelled along the crystal c≡z axis with a set of home-built superconducting magnet coils.

The various optical transitions of a single ¹⁷¹Yb qubit are employed for state readout and initialisation (FIG. 5 a ). Optical addressing of the A transition for readout is established with a continuous-wave (CW) titanium sapphire (Ti:Sapph) laser (M² Solstis) which is frequency-stabilised to a high-finesse reference cavity (Stable Laser Systems) using Pound-Drever-Hall locking [43]. The laser double-passes through two freespace acousto-optic-modulator (AOM) setups leading to single-photon level extinction of the input beam, and pulse generation with ≈10 ns rise times. A second CW external cavity diode laser (Toptica DL-Pro) is used to address the F transition during initialisation. The laser passes through an identical AOM setup and is frequency stabilised via offset-frequency locking to the Ti:Sapph.

The light output from the cavity is separated from the input with a 99:1 fibre beamsplitter, and passed through a single AOM which provides time-resolved gating of the light to prevent reflected laser pulses from saturating the detector. The light is then sent to a tungsten-silicide superconducting nanowire single photon detector (SNSPD) (Photonspot) which also sits on the still-plate of the cryostat. Photon detection events are subsequently timetagged and histogrammed (Swabian Timetagger 20).

Microwave pulses to control the ground-state qubit transition (6.75 MHz) and square-wave RF to generate the ¹⁷¹Yb⁵¹V interaction (100-300 kHz) are directly synthesised with an arbitrary waveform generator (Tektronix 5204AWG) and amplified (Amplifier Research 10U1000). A second microwave path is used for the excited state microwave control (3.4 GHz) necessary for qubit initialisation. The control pulses are generated by switching the output of a signal generator (SRS SG386) and amplifying (Minicircuits ZHL-16 W-43-S+). The two microwave signal paths are combined with a diplexer (Marki DPXN2) and sent into the fridge to the device. A gold coplanar waveguide fabricated on the YVO₄ surface enables microwave driving of the ions.

b. ¹⁷¹Yb Initialisation, Readout and Experiment Sequence

At the 500 mK experiment operating temperature and at zero magnetic field, the equilibrium ¹⁷¹Yb population is distributed between the |aux_(g)

, 0_(g)

and |1_(g)

states (FIG. 5 a ) All experiments start by initialising the single ¹⁷¹Yb, ion into |0_(g)

via a two-stage protocol [15]. Firstly the |aux_(g)

state is emptied with a series of 3 μs pulses applied to the optical F transition each followed by a 3 μs wait period. When the ¹⁷¹Yb ion is successfully excited from |aux_(g)

to |1_(g)

, the population in |1_(g)

will preferentially decay to |0_(g)

during the wait time via the cavity-enhanced E transition. Subsequently, the |1_(g)

state is also emptied by applying an optical π pulse to the A transition followed by a microwave π pulse to the f_(g) transition in rapid succession, which similarly leads to excitation from |1_(g)

to |1_(e)

and decay into |0_(g)

. This process is repeated several times for improved fidelity.

Readout of the ¹⁷¹Yb|1_(g)

state is performed by applying a series of 100π pulses to the A transition, each of which is followed by a 10 μs photon detection window. This process is enabled by the cyclic nature of the A transition. To read out the |0_(g)

population we apply an additional π pulse to swap the |0_(g)

↔|1_(g)

populations before performing the same optical readout procedure.

FIG. 5 c shows an exemplary pulse sequence used to store and retrieve a superposition state from the register consisting of four ⁵¹V lattice ions. The sequence starts with initialisation of the ¹⁷¹Yb qubit into |0_(g)

and the ⁵¹V register into |0_(v)

=|±7/2

^(⊗4). A series of ZenPol polarisation operations are interleaved with ¹⁷¹Yb re-initialisation sequences and alternate between ω_(b) and ω_(c) transition control to sequentially polarize the spin-7/2 ⁵¹V register towards the |±7/2

level. After the initialization sequence, a single π/2 pulse is applied to the ¹⁷¹Yb qubit to prepare a superposition state. Subsequently, the state is transferred to the ⁵¹V register using a swap operation resonant with the ω_(c) transition as detailed in the main text. After a variable wait time, the superposition state is retrieved with a second swap gate and measured in the x-basis via π/2 pulse followed by optical readout on the A transition as detailed above.

c. ZenPol Sequence

Consider a system of a single ¹⁷¹Yb qubit coupled to four neighbouring nuclear spin-7/2 ⁵¹V ions. This hybrid spin system is described by the effective Hamiltonian (setting ℏ=1):

$\hat{H} = {{{\Delta(t)}{\hat{\overset{\sim}{S}}}_{z}} + {\sum\limits_{i \in {register}}{Q\left( {\hat{I}}_{z}^{(i)} \right)}^{2}} + {\sum\limits_{i \in {register}}{{{\hat{\overset{\sim}{S}}}_{z}\left\lbrack {B_{z}^{OH} + {B^{RF}(t)}} \right\rbrack}\left\lbrack {{a_{x}{\hat{I}}_{x}^{(i)}} + {a_{z}{\hat{I}}_{z}^{(i)}}} \right\rbrack}}}$

where Δ(t)=γ_(z) ²(B_(π) ^(OH)+B^(RF)(t))²/2ω₀₁ is the effective energy shift due to both z-directed nuclear Overhauser (B_(π) ^(OH)) and external RF(B^(RF)(t)) magnetic fields, ω₀₁/2π=675 MHz is the ¹⁷¹Yb qubit transition frequency, γ_(z)/27=8.5 MHz/G is the ¹⁷¹Yb ground-state longitudinal gyromagnetic ratio, Q/2π=165 kHz is the ⁵¹V register nuclear quadrupole splitting, {tilde over (Ŝ)}_(x) is the ¹⁷¹Yb qubit operator along the z-axis, Î_(x,z) are the ⁵¹V spin-7/2 operators along the x- and z-axis, and a_(x,z) are the effective coupling strengths between ¹⁷¹Yb and ⁵¹V along the x- and z-axes. See Supplementary Information for a detailed derivation of this effective Hamiltonian.

As discussed in the main text, polarisation of the ⁵¹V register and preparation of collective spin-wave states relies on induced polarisation transfer from the ¹⁷¹Yb to ⁵¹V and is achieved via periodic driving of the ¹⁷¹Yb qubit. Specifically, periodic pulsed control can dynamically engineer the original Hamiltonian (equation (7)) to realize effective spin-exchange interaction between ¹⁵¹Yb and ⁵¹V ions of the form, {tilde over (Ŝ)}₊Î⁻+{tilde over (Ŝ)}⁻Î₊, in the average Hamiltonian picture [36],[44] One example of such a protocol is the recently developed PulsePol sequence [35], however, it relies on states with a constant, nonzero magnetic dipole moment and therefore cannot be used in our system since the ¹⁷¹Yb qubit has no intrinsic magnetic dipole moment. Motivated by this approach, we have developed a variant of the PulsePol sequence that accompanies a square-wave RF field synchronized with the sequence (FIG. 8 a ). The base sequence has a total of 8 free-evolution intervals with equal duration (τ/4) defined by periodically spaced short pulses and is repeatedly applied to ¹⁷¹Yb. Following the sequence design framework presented in Ref. [36], we judiciously choose the phase and ordering of the constituent π/2 and π pulses such that the resulting effective interaction has spin-exchange form with strength proportional to the RF magnetic field amplitude (B^(RF)), whilst decoupling from interactions induced by the Overhauser field (B_(π) ^(OH)). The sequence was designed to cancel detuning induced by both of these fields and to retain robustness against pulse rotation errors to leading order. We term this new sequence ‘ZenPol’ for ‘zero first-order Zeeman nuclear-spin polarisation’.

To understand how the ZenPol sequence works, one can consider a toggling-frame transformation of the ¹⁷¹Yb spin operator along the quantisation axis ({tilde over (Ŝ)}_(z,tog)(t)): we keep track of how this operator is transformed after each preceding pulse. For example, the first π/2 pulse around the y-axis transforms {tilde over (Ŝ)}_(π) into −{tilde over (Ŝ)}_(x) and the subsequent π pulse around the y-axis transforms −{tilde over (Ŝ)}_(x) into +{tilde over (Ŝ)}_(x). Over one sequence period, the toggling frame transformation generates a time-dependent Hamiltonian Ĥ_(tog)(t) that is piecewise constant for each of 8 free-evolution intervals, which can be expressed as

${{\hat{H}}_{tog}(t)} = {{{\Delta(t)}\left\lbrack {{{f_{x}^{OH}(t)}{\hat{\overset{\sim}{S}}}_{x}} + {{f_{y}^{OH}(t)}{\hat{\overset{\sim}{S}}}_{y}}} \right\rbrack} + {\sum\limits_{i \in {register}}{Q\left( {\hat{I}}_{z}^{(i)} \right)}^{2}} + {\sum\limits_{i \in {register}}{{B_{z}^{OH}\left\lbrack {{{f_{x}^{OH}(t)}{\hat{\overset{\sim}{S}}}_{x}} + {{f_{y}^{OH}(t)}{\hat{\overset{\sim}{S}}}_{y}}} \right\rbrack}\left\lbrack {{a_{x}{\hat{I}}_{x}^{(i)}} + {a_{z}{\hat{I}}_{z}^{(i)}}} \right\rbrack}} + {\sum\limits_{i \in {register}}{{B^{RF}\left\lbrack {{{f_{x}^{RF}(t)}{\hat{\overset{\sim}{S}}}_{x}} + {{f_{y}^{RF}(t)}{\hat{\overset{\sim}{S}}}_{y}}} \right\rbrack}\left\lbrack {{a_{x}{\hat{I}}_{x}^{(i)}} + {a_{z}{\hat{I}}_{z}^{(i)}}} \right\rbrack}}}$

Here f_(x,y) ^(OH)(t) describes the time-dependent modulation of the ¹⁷¹Yb z-spin operator ({tilde over (Ŝ)}_(z.tog)(t)=f_(x) ^(OH)(t){tilde over (Ŝ)}_(x)+f_(y) ^(OH)(t){tilde over (Ŝ)}_(y)) (FIG. 8 a ). Note that f_(π) ^(OH)(t)=0 for all intervals. Since the externally-applied squarewave RF field is constant for each half-sequence period, we can replace B^(RF)(t) with the amplitude B^(RF) and transfer the time dependence to f_(x,y) ^(OH) by applying sign flips, thus leading to redefined modulation functions f_(x,y) ^(RF) (FIG. 8 a ).

The spin-7/2 ⁵¹V ion exhibits three distinct transitions at frequencies ω_(a,b,c) (FIG. 1 b ) In the following, we consider an effective spin-1/2 system for the ⁵¹V ions using the ω_(c) manifold {|Î

=|=5/2

, |↓

=|±7/2

}, with

${{\hat{\overset{\sim}{I}}}_{x} = {\frac{1}{2}\left( ❘\uparrow \right)\left( {\left. \downarrow ❘ \right. + \left. ❘\downarrow \right.} \right)\left( \uparrow ❘ \right)}},{{\hat{\overset{\sim}{I}}}_{y} = {{\frac{1}{2i}\left( ❘\uparrow \right)\left( {\left. \downarrow ❘ \right. - \left. ❘\downarrow \right.} \right)\left( \uparrow ❘ \right){and}{\hat{\overset{\sim}{I}}}_{z}} = {\frac{1}{2}\left( ❘\uparrow \right)\left( \uparrow{❘ - ❘}\downarrow \right){\left( \downarrow ❘ \right).}}}}$

In a rotating frame with respect to the target frequency ω_(c), the nuclear spin operators become {tilde over (Î)}_(x)→{tilde over (Î)}_(x) cos(ω_(c)t)+{tilde over (Î)}_(y) sin(ω_(c)t) and {tilde over (Î)}π→{tilde over (Î)}_(x). Thus, the leading- order average Hamiltonian,

Ĥ_(avg) = ?dtĤ_(tog)(t), ?indicates text missing or illegible when filed

in the rotating frame is given by:

${\hat{H}}_{avg} = {\sum\limits_{i \in {register}}{\frac{a_{x}\sqrt{7}}{2\tau}{\int_{0}^{2\tau}{{dt}{\left\{ {{{B_{z}^{OH}\left\lbrack {{{f_{x}^{OH}(t)}{\hat{\overset{\sim}{S}}}_{x}} + {{f_{y}^{OH}(t)}{\hat{\overset{\sim}{S}}}_{y}}} \right\rbrack}\left\lbrack {{{\hat{\overset{\sim}{I}}}_{x}^{(i)}{\cos\left( \text{?} \right)}} + {{\hat{\overset{\sim}{I}}}_{y}^{(i)}{\sin\left( \text{?} \right)}}} \right\rbrack} + {{B^{RF}\left\lbrack {{{f_{x}^{RF}(t)}{\hat{\overset{\sim}{S}}}_{x}} + {{f_{y}^{RF}(t)}{\hat{\overset{\sim}{S}}}_{y}}} \right\rbrack}\left\lbrack {{{\hat{\overset{\sim}{I}}}_{x}^{(i)}{\cos\left( \text{?} \right)}} + {{\hat{\overset{\sim}{I}}}_{y}^{(i)}{\sin\left( \text{?} \right)}}} \right\rbrack}} \right\}.}}}}}$ ?indicates text missing or illegible when filed

Here, various terms are excluded as they time average to zero (rotating-wave approximation). The √{square root over (7)} prefactor comes from mapping the original spin-7/2 operators to the effective spin-1/2 ones. Additionally, the energy shift induced by B_(π) ^(TH) and time-dependent B^(RF) is cancelled since we are using square-wave RF. The Fourier transforms of the modulation functions f_(x,y)(t), termed the filter functions [45], directly reveal resonance frequencies at which equation (9) yields non-zero contributions (FIG. 8 b ). Resonant interactions with strength proportional to the nuclear Overhauser field are achieved at sequence periods 2τ which satisfy

${\frac{1}{2\tau} = \frac{\text{?}}{2{\pi \times 2}}},\frac{\text{?}}{2{\pi \times 4}},\frac{\text{?}}{2{\pi \times 6}},{\cdots;}$ ?indicates text missing or illegible when filed

interactions proportional to the RF field occur at sequence periods satisfying

${\frac{1}{2\tau} = \frac{\text{?}}{2{\pi \times 1}}},\frac{\text{?}}{2{\pi \times 3}},\frac{\text{?}}{2{\pi \times 5}},{\cdots.}$ ?indicates text missing or illegible when filed

These two sets of resonances occur at different values of 2τ, hence we can preferentially utilise the coherent, RF-induced interactions whilst decoupling from those induced by the randomised Overhauser field. This is experimentally demonstrated in FIG. 2 b where the RF-induced resonances are spectrally resolved. In this measurement the linewidth of the register resonances are limited by that of the filter function. We also note that the transition cannot be independently addressed by the ZenPol sequence due to the multiplicity of the three ⁵¹V transitions determined by the quadratic Hamiltonian (ω_(a)=ωb/2=ω_(c)/3).

We use the RF-driven resonance identified at

$\frac{1}{\text{?}} = \frac{\text{?}}{2{\pi \times 5}}$ ?indicates text missing or illegible when filed

by setting the free-evolution interval to

$\frac{\text{?}}{\text{?}} = {\frac{\text{?}}{\text{?}}.}$ ?indicates text missing or illegible when filed

Under this resonance condition, the average Hamiltonian (equation (9)) is simplified to

${\hat{H}}_{avg} = {{{- \sqrt{7}}\left( \frac{1 + \sqrt{2}}{5\pi} \right)a_{x}{B^{RF} \times {\sum\limits_{i \in {register}}\left( {{\left( {{\hat{\overset{\sim}{S}}}_{x} + {\hat{\overset{\sim}{S}}}_{y}} \right){\hat{\overset{\sim}{I}}}_{x}^{(i)}} + {\left( {{- {\hat{\overset{\sim}{S}}}_{x}} + {\hat{\overset{\sim}{S}}}_{y}} \right){\hat{\overset{\sim}{I}}}_{y}^{(i)}}} \right)}}} = {{{- \sqrt{7}}\left( \frac{\sqrt{2} + 2}{5\pi} \right)a_{x}B^{RF}{\sum\limits_{i \in {register}}\left( {{{\hat{\overset{\sim}{S}}}_{x}^{\prime}{\hat{\overset{\sim}{I}}}_{x}^{(i)}} + {{\hat{\overset{\sim}{S}}}_{y}^{\prime}{\hat{\overset{\sim}{I}}}_{y}^{(i)}}} \right)}} = {\text{?}B^{RF}{\sum\limits_{i \in {register}}{\left( {{{\hat{\overset{\sim}{S}}}_{+}^{\prime}{\hat{\overset{\sim}{I}}}_{-}^{(i)}} + {{\hat{\overset{\sim}{S}}}_{-}^{\prime}{\hat{\overset{\sim}{I}}}_{+}^{(i)}}} \right).}}}}}$ ?indicates text missing or illegible when filed

Here, going from the first to the second line, we change the local ¹⁷¹Yb basis by rotating 45 degrees around the z axis such that {tilde over (Ŝ)}_(x)′=({tilde over (Ŝ)}_(x)+{tilde over (Ŝ)}_(y))/√{square root over (2)}, {tilde over (Ŝ)}_(y)′=(−{tilde over (Ŝ)}_(x)+{tilde over (Ŝ)}_(y))/√{square root over (2)}, and from the second to the third line, {tilde over (Ŝ)}₊′={tilde over (Ŝ)}_(x)′±i{tilde over (Ŝ)}_(y)′ and {tilde over (Î)}₊′={tilde over (Î)}_(x)′±i{tilde over (Ŝ)}_(y) are used. We define the coefficient b_((k,ω) _(j) ₎ which determines the interaction strength for the k^(th) resonance addressing transition ω_(i) (for example, b_((5,ω) _(c) ₎=−√{square root over (7)}(√{square root over (2)}+2)a_(x)/10π). In the discussion in the first example above, we omitted the primes on the ¹⁷¹Yb spin operators for the sake of notational simplicity. The same analysis can be performed for other transitions, yielding a similar spin exchange Hamiltonian, albeit with different interaction 9) strength.

d. Direct Drive Gates for ⁵¹V Register

Performing dynamical decoupling on the register requires selective driving of the froze-core ⁵¹V nuclear spins without perturbing the bath and is achieved through a two-fold mechanism. Firstly we initialise the ¹⁷¹Yb qubit into |0_(g)

and apply a sinusoidal z-directed RF magnetic field at ω_(c)/2π=991 kHz through the coplanar waveguide to induce an oscillating ¹⁷¹Yb magnetic dipole moment (FIG. 11 a ) This generates an x-directed field component at each ⁵¹V spin, where the driving Hamiltonian is given by Ĥ_(drive)=μ_(N)q_(vx)A_(x)B_(π) ^(osc) sin(ω_(c)t)Ĩ_(x) with A_(x)=3 ln μ₀γ_(z) ²/8πr³ω₀₁. Here, μ_(N) is the nuclear magneton, g_(vx) is the ⁵¹V directed g-factor, B_(z) ^(osc) is the sinusoidal RF magnetic field amplitude, Î_(x) is the nuclear spin-7/2 operator along the x-axis, {l,n} are the {x,z} directional cosines of the ¹⁷¹Yb-⁵¹V displacement vector, μ₀ is the vacuum permittivity, and r is the ¹⁷¹Yb-⁵¹V ion distance (Supplementary Information). The lattice symmetry of the host leads to equidistant spacing of the four proximal ⁵¹V spins from the central ¹⁷¹Yb qubit allowing homogeneous coherent driving of all register spins.

In this direct driving scheme, we note that the effect of B_(z) ^(osc) is amplified by a factor of A_(x)≈6.7 for the frozen core register spins at a distance of r=3.9 Å (Supplementary Information). Crucially, the amplification factor scales as A_(x)∝1/r³ with distance r from the ¹⁷¹Yb qubit, leading to a reduced driving strength for distant ⁵¹V bath spins. Moreover, the transition frequency of the bath, ω_(c) ^(bath)/2π=1028 kHz, is detuned by 37 kHz from that of the register, ω_(c)/2π=991 kHz, further weakening the bath interaction due to off-resonant driving provided that the Rabi frequency is less than the detuning.

In a rotating frame at frequency ω_(c), the driving Hamiltonian Ĥ_(drive) gives rise to Rabi oscillation dynamics of the register spins within the ω_(c) manifold, {|↑

=|±5/2

, |↓

=|±7/2

}. To calibrate ⁵¹V π pulse times, we initialise the register into |0_(v)

=|↓↓↓↓

, drive the register for variable time, and read out the |0_(v)

population by preparing the ¹⁷¹Yb qubit in |1_(g)

and applying a swap gate to the ω_(c) transition. If the final ⁵¹V spin state is in |↓

(|↑

) the swap will be successful (unsuccessful) and the ¹⁷¹Yb qubit will end up in |0_(g)

(|1_(g)

). Using this method, we induce resonant Rabi oscillations of the register at a Rabi frequency of 2π×(7.65±0.05) kHz (blue markers, FIG. 11 c ) which exhibit exponential decay on a 280±30 μs timescale, limited by dephasing caused by the fluctuating ¹⁷¹Yb Knight field. This can be decoupled using motional narrowing techniques whereby we periodically apply π pulses to the ¹⁷¹Yb every 6 μs during the drive period. In order to drive the ⁵¹V spins in a phase-continuous manner, we compensate for the inversion of the ¹⁷¹Yb magnetic dipole moment after each π pulse by applying a π phase shift to the sinusoidal driving field (FIG. 11 b ). This leads to an extended 1/e Gaussian decay time of 1040±70 μs (red markers, FIG. 11 c ). The arrow in FIG. 11 c indicates the 69 μs ⁵¹Vπ pulse time used for dynamical decoupling. In contrast to the spin-preserving exchange interaction, this direct drive protocol provides independent, local control of the four ¹⁷¹V spins with no constraints on the number of excitations, thereby coupling the ⁵¹V register to states outside the two-level manifold spanned by |0_(v)

and |W_(v)

. For example, at odd multiple π times, we find

$\left. \left. {\left. {\left. \left. {❘0_{v}} \right\rangle\rightarrow \right.❘}\uparrow\uparrow\uparrow\uparrow \right\rangle{❘W_{v}}} \right\rangle\rightarrow\frac{\left. \left. {\left. {\left. {\left. \left( {❘\left. \downarrow\uparrow\uparrow\uparrow \right.} \right. \right\rangle + {❘\left. \uparrow\downarrow\uparrow\uparrow \right.}} \right\rangle + {❘\left. \uparrow\uparrow\downarrow\uparrow \right.}} \right\rangle + {❘\left. \uparrow\uparrow\uparrow\downarrow \right.}} \right\rangle \right)}{2} \right.,$

both of which contain more than a single excitation. For this reason, we use an even number of π pulses in our decoupling sequences to always return the ⁵¹V register to the memory manifold prior to state retrieval.

e. Population Basis Measurements

We developed a sequential tomography protocol [39] to read out the populations of the joint ¹⁷¹Yb-⁵¹V density matrix ρ in the effective four-state basis, {|0_(g)0_(v)

, |0_(g)W_(v)

, |1_(g)0_(v)

, |1_(g)W_(v)

}. This is achieved using two separate sequences: Readout sequence 1 and Readout sequence 2, applied alternately, which measure the {|0_(g)0_(v)

, |0_(g)W_(v)

} and {|1_(g)0_(v)

, |1_(g)W_(v)

} populations respectively. As shown in FIG. 13 a , these sequences are distinguished by the presence (absence) of a single π pulse applied to the ¹⁷¹Yb qubit at the start of the sequence. This is followed by a single optical readout cycle on the A transition; results are post-selected on detection of a single optical photon during this period. Hence the presence (absence) of the first π pulse results in |0_(g)

(|1_(g)

) state readout after post selection. Furthermore, in all post-selected cases the ¹⁷¹Yb qubit is initialised to |1_(g)

by taking into account this conditional measurement outcome. Subsequently, an unconditional π pulse is applied to the ¹⁷¹Yb, preparing it in |0_(g)

and a swap gate is applied, thereby transferring the ⁵¹V state to the ¹⁷¹Yb. Finally, we perform single-shot readout of the ¹⁷¹Yb state according to the protocol developed in [15]. Specifically, we apply two sets of 100 readout cycles to the A transition separated by a single π pulse which inverts the ¹⁷¹Yb qubit population. The ⁵¹V state is ascribed to |W_(v)

(|0_(v)

) if ≥1(0) photons are detected in the second readout period and 0(≥1) photons are detected in the third. The possible photon detection events and state attributions are summarized in FIG. 13 b.

This protocol was demonstrated by characterizing the state preparation fidelities of the four basis states, the measured histograms are presented in FIG. 13 c alongside the respective gate sequences used for state preparation. The resulting uncorrected (corrected) preparation fidelities for these four basis states are:

_(|0) _(g) ₀ _(v)

=0.79±0.01(0.82±0.02),

_(|0) _(g) _(W) _(v)

=0.50±0.02(0.64±0.02),

_(|1) _(g) ₀ _(v)

=0.79±0.01(0.82±0.02),

_(|1) _(g) _(W) _(v)

=0.50±0.02(0.61±0.02)

The reduced fidelity of |0_(g)W_(v)

and |1_(g)W_(v)

relative to |0_(g)0_(v)

and |1_(g)0_(v)

arises from the swap gate used for the |W_(v)

state preparation. Finally, we characterized the fidelity of the maximally entangled ¹⁷¹Yb-⁵¹V Bell state,

${\left. ❘\Psi^{+} \right\rangle = {\left. \frac{1}{\sqrt{2}}\text{(}1_{g}\text{?} \right\rangle - {\left. i❘0_{g}W_{v} \right\rangle)}}},$ ?indicates text missing or illegible when filed

prepared using a single √{square root over (swap)} gate as described in the first example (FIG. 13 d ). The corresponding uncorrected (corrected) populations for the four basis states, p_(ij)(c_(ij)) denoted are:

p ₀₀≡

0_(g)0_(v)|ρ|0_(g)0_(v)

=0.16±0.01(c ₀₀=0.07±0.02),

p ₀₁≡

0_(g) W _(v)|ρ|0_(g) W _(v)

=0.326±0.01(c ₀₁=0.41±0.02),

p ₀₁≡

1_(g)0_(v)|ρ|1_(g)0_(v)

=0.40±0.02(c ₁₀=0.41±0.02),

p ₁₁≡

1_(g) W _(v)|ρ|1_(g) W _(v)

=0.12±0.01(c ₁₁=0.11±0.01),

f. Swap Gate Fidelity Correction

Since ¹⁷¹Yb readout fidelity is >95% [15], the dominant error introduced during the population basis measurements arises from the swap gate. Its fidelity in the population basis was measured by preparing either the |0_(g)0_(v)

state (zero spin excitations) or the |1_(g)0_(v)

state (single spin excitation) and applying two consecutive swap gates such that the system is returned to the initial state. By comparing the ¹⁷¹Yb population before (p_(pre)) and after (p_(post)) the two gates are applied, fidelity estimates can be extracted independently from the ⁵¹V state initialisation. Assuming the swap and swapback processes are symmetric, a gate fidelity

_(sw)=√{square root over ((1−2p_(post))/(1−2p_(pre)))} is obtained. This quantity is measured for zero spin excitations leading to

_(sw,0)=0.83 and with a single spin excitation leading to

_(sw,1)=0.52.

When measuring the joint ¹⁷¹Yb-⁵¹V populations {p₀₀, p₀₁, p₁₀, p₁₁}, these fidelities can be used to extract a set of corrected populations {c₀₀, c₀₁, c₁₀, c₁₁} according to the method described in [46,47] using

$\begin{pmatrix} c_{11} \\ c_{10} \\ c_{01} \\ c_{00} \end{pmatrix} = {E^{- 1}\begin{pmatrix} p_{11} \\ p_{10} \\ p_{01} \\ p_{00} \end{pmatrix}}$

where

$E = {\frac{1}{2}{\begin{pmatrix} {1 + \mathcal{F}_{{sw},1}} & {1 - \mathcal{F}_{{sw},0}} & 0 & 0 \\ {1 - \mathcal{F}_{{sw},1}} & {1 + \mathcal{F}_{{sw},0}} & 0 & 0 \\ 0 & 0 & {1 + \mathcal{F}_{{sw},1}} & {1 - \mathcal{F}_{{sw},0}} \\ 0 & 0 & {1 - \mathcal{F}_{{sw},1}} & {1 + \mathcal{F}_{{sw},0}} \end{pmatrix}.}}$

A similar approach to correct the √{square root over (swap)} gate was used to read out the Bell state coherence.

6. Supplementary Example Derivations for Interactions and Hamiltonians Described Herein

a. ¹⁷¹Yb-⁵¹-V Interactions

(i). Ground State ¹⁷¹Yb Hamiltonian

The effective spin-1/2 Hamiltonian for the ^(z)F_(7/2)(0)¹⁷¹Yb³⁺ ground state is given by [1]:

Ĥ _(eff) =μBB·g·Ŝ+Î _(Vh) ·A·Ŝ

where B is the magnetic field, Ŝ and Î_(Yb) are vectors of ¹⁷¹Yb electron and nuclear spin-1/2 operators respectively and we neglect the nuclear Zeeman term. The uniaxial ground state g tensor is given by:

$g = {\begin{pmatrix} g_{x} & 0 & 0 \\ 0 & g_{x} & 0 \\ 0 & 0 & g_{z} \end{pmatrix} = \begin{pmatrix} 0.85 & 0 & 0 \\ 0 & 0.85 & 0 \\ 0 & 0 & {- 6.08} \end{pmatrix}}$

which is a uniaxial tensor with the extraordinary axis parallel to the c-axis of the crystal and the two ordinary axes aligned with the crystal a-axes. The ground state A tensor is given by:

$A = {2{\pi \times \begin{pmatrix} 0.675 & 0 & 0 \\ 0 & 0.675 & 0 \\ 0 & 0 & {- 4.82} \end{pmatrix}}{GHz}}$

FIG. 5 a shows the zero magnetic field energy level structure with hybridised ¹⁷¹Yb electron-nuclear spin eigenstates. Note that the zero-field ¹⁷¹Yb qubit states, |0_(g)

and |1_(g)

, have no magnetic dipole moment. See [1] for more details. Throughout this work we adopt ℏ=1 convention.

(ii). Local Nuclear Spin Environment

The ¹⁷¹Yb³⁺ ion substitutes for yttrium in a single site of the YVO₄ crystal, furthermore naturally abundant Y and V contain 99.8% ⁵¹V and 100% ⁸⁹Y isotopes. Hence each ¹⁷¹Yb ion experiences a near-identical nuclear spin environment. The ⁵¹V ions have nuclear spin-7/2 leading to electric quadrupole interactions that cause a zero-field splitting. The resulting zero-field energy level structure of the bath is given by:

Ĥ _(v) =Q _(bath) Î _(x) ²

with Q_(bath)/2π=171 kHz measured using nuclear magnetic resonance (NMR) on bulk YVO₄ crystals [2] and Î_(x) the ⁵¹V nuclear spin-7/2 spin operator along the c≡z axis. Note that the local ⁵¹V register ions surrounding the ¹⁷¹Yb qubit experience a frozen-core detuning as discussed in the main text, leading to a smaller quadrupolar splitting with Q/2π=165 kHz. The energy level structure of these register ions is shown in FIG. 1 b . The ⁸⁹Y ion, on the other hand, has no zero-field structure.

The positions of the six nearest ⁵¹V ions are tabulated below, where r=[xyz] is the ¹⁷¹Yb-⁵¹V position vector with magnitude r and direction cosines {l,m,n}.

r x y z ⁵¹V ion # Shell (Å) (Å) (Å) (Å) l m n 1 1^(st) 3.1 0 0 −3.1 0 0 −1 2 1^(st) 3.1 0 0 3.1 0 0 1 3 2^(nd) 3.9 0 −3.6 1.6 0 −0.91 0.40 4 2^(nd) 3.9 0 3.6 1.6 0 0.91 0.40 5 2^(nd) 3.9 −3.6 0 −1.6 −0.91 0 −0.40 6 2^(nd) 3.9 3.6 0 −1.6 0.91 9 −0.40

Note that the two nearest ⁵¹V ions (1 and 2) are located directly above and below the ¹⁷¹Yb qubit along the z-axis, due to their positions they cannot be driven by the induced ¹⁷¹Yb magnetic dipole moment and thus belong to the bath (Supplementary Information Section I C). In contrast, ions 3-6 are symmetrically positioned in the lattice with non-zero x/y and z coordinates, forming the frozen-core register spins utilized as a quantum memory. The ⁵¹V ions have a uniaxial g-tensor with form [3]:

$g_{v} = \begin{pmatrix} g_{vx} & 0 & 0 \\ 0 & g_{vx} & 0 \\ 0 & 0 & g_{vz} \end{pmatrix}$

(iii) ¹⁷¹Yb-⁵¹V Interactions

The magnetic dipole-dipole interaction between the ¹⁷¹Yb qubit and a single ⁵¹V ion can be described by the following Hamiltonian:

${\hat{H}}_{dd} = {\frac{\mu_{0}}{4\pi}\left\lbrack {\frac{\mu_{Yb} \cdot \mu_{V}}{r^{3}} - \frac{3\left( {\mu_{Yb} \cdot r} \right)\left( {\mu_{V} \cdot r} \right)}{r^{5}}} \right\rbrack}$

where μ_(Vh)=−μ_(B)g·Ŝ, μ_(V)=μ_(N)g_(V)·Î (note that Î is a vector of ⁵¹V nuclear spin operators), μ_(B) is the Bohr magneton, μ_(N) is the nuclear magneton, μ₀ is the vacuum permeability and r is the ¹⁷¹Yb ⁵¹V displacement vector with magnitude r. Due to the highly off-resonant nature of the ¹⁷¹Yb⁵¹V interaction, a secular approximation would be appropriate. To first order, however, all secular terms involving the ¹⁷¹Yb qubit basis are zero, i.e.,

0_(g)|Ĥ_(dd)|0_(g)

=0,

1_(g)|Ĥ_(dd)|1_(g)

=0

To proceed, consider that second-order effects which generally scale as ˜g²/ΔE, where ΔE is the energy separation between a pair of unperturbed eigenstates. By taking into account the fact that g_(z) is roughly 7 times larger than g_(x),g_(y) and Ŝ_(π) terms in Ĥ_(dd) mix |0_(g)

and |1_(g)

with small ΔE whereas, Ŝ_(x) and Ŝ_(y) mix the ¹⁷¹Yb qubit states and |aux g

with large ΔE, we restrict our consideration to the Ŝ_(π) terms in Ĥ_(dd):

${\hat{H}}_{dd} \approx {\frac{\mu_{0}\mu_{B}\mu_{N}g_{z}}{4\pi r^{3}}{{\hat{S}}_{z}\left\lbrack {{3\ln g_{vx}{\hat{I}}_{x}} + {3{mng}_{vx}{\hat{I}}_{y}} + {\left( {{3n^{2}} - 1} \right)g_{vz}{\hat{I}}_{z}}} \right\rbrack}}$

where {l,m,n} are direction cosines of the ¹⁷¹Yb⁵¹V displacement vector. Note that the Ŝ_(x) operator is the electron spin-1/2 operator defined as Ŝ_(z)=½(|0_(g)

1_(g)|+|1_(g)

0_(g)|) in the basis of the hybridised eigenstates of the ¹⁷¹Yb qubit.

(iv) Nuclear Overhauser Field

As discussed in the first example, the ⁵¹V spins can be divided into two ensembles: register spins and bath spins. The bath spins comprise ⁵¹V ions which are not driven by the ¹⁷¹Yb qubit for the following two reasons:

1 Ions which aren't driven due to position: certain ions (such as 1 and 2 in the above table) only interact via an Ising-type Ŝ_(x)Î_(x) Hamiltonian. Hence the ¹⁷¹Yb qubit cannot be used to drive transitions between the ⁵¹V z-quantised quadrupole levels.

2 Ions which aren't driven due to detuning: As observed in the ZenPol spectra (FIG. 2 b in the main text), more distant spins are spectrally separated from the nearby ions comprising the register.

It is assumed that the bath spins are in an infinite-temperature mixed state: ρ_(V)=1_(V)/Tr{1_(V)}, where 1_(V) is the identity matrix in the Hilbert space for the bath spins. In the mean field picture, their effect on the ¹⁷¹Yb can be approximated as a classical fluctuating magnetic field, commonly termed the nuclear Overhauser field. As mentioned previously, since g_(z) ²>>g_(x,y) ², the z-component of the Overhauser field is dominant, given by

${B_{z}^{OH} = {\sum\limits_{i \in {bath}}{\frac{\mu_{0}\mu_{N}g_{vz}}{4{\pi\left( r^{(i)} \right)}^{3}}\left( {{3\left( n^{(i)} \right)^{2}} - 1} \right)m_{I}^{(i)}}}},$

where r_((i)) and n^((i)) are the distance and z-direction cosine between the ¹⁷¹Yb and i^(th) bath spin, and m_(l) ^((i))∈{−7/2, −5/2, . . . , 5/2, 7/2} is the nuclear spin projection at site i. Note that B_(π) ^(OH) is randomly fluctuating due to the stochastic occupation of the 8 possible |m_(I)

states, however, it is quasi-static on the timescale of our control sequences, hence we do not label the time dependence. The nuclear Overhauser field generates some weak mixing between |0_(g)

and |1_(g)

leading to perturbed eigenstates |{tilde over (0)}_(g)

and |{tilde over (1)}_(g)

which have a small, induced, z-directed dipole moment. x These states have the form

${\left. ❘{\overset{\sim}{0}}_{g} \right\rangle = {\left. ❘0_{g} \right\rangle + \left. \frac{\gamma_{z}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}{2\omega_{01}}❘1_{g} \right\rangle}}{\left. ❘{\overset{\sim}{1}}_{g} \right\rangle = {\left. ❘1_{g} \right\rangle - \left. \frac{\gamma_{z}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}{2\omega_{01}}❘0_{g} \right\rangle}}$

where γ_(z)=g_(z)μ_(B) is the longitudinal gyromagnetic ratio of the ¹⁷¹Yb qubit and ω₀₁/2π=675 MHz is the unperturbed ¹⁷¹Yb|0_(g)

↔|1_(g)

transition frequency. Here we have added the effect of an externally applied, z-directed, square wave RF magnetic field B^(RF)(t) with amplitude B^(RF) used in the ZenPol sequence; note that this field is piecewise constant for each half-sequence period, hence the time dependence corresponds to periodic flips between ±B^(RF). In addition, these fields induce a detuning of the ¹⁷¹Yb|0_(g)

↔|1_(g)

transition, which can be calculated using second-order perturbation theory as Δ(t)=γ_(x) ²(B_(π) ^(OH)+B^(RF)(t))²/2ω₀₁.

(v) Interaction with Register Ions

We postulate that the second nearest shell of four ⁵¹V ions (ions 3-6 in the table above) comprise the register. These four ions are equidistant from the ¹⁷¹Yb and interact via both an Ŝ_(x)Î_(x) term and Ŝ_(z)Î_(x) or Ŝ_(z)Î_(y) terms. To identify an effective interaction Hamiltonian in the perturbed basis {|{tilde over (0)}_(g)

, |{tilde over (1)}_(g)

}, only secular matrix elements of Ĥ_(dd) (equation (S7)) are considered:

${\hat{\overset{\sim}{H}}}_{dd} = {{\left\langle {\overset{\sim}{0}}_{g}❘{\hat{H}}_{dd}❘{\overset{\sim}{0}}_{g} \right\rangle\left. ❘{\overset{\sim}{0}}_{g} \right\rangle\left\langle {\overset{\sim}{0}}_{g}❘ \right.} + {\left\langle {\overset{\sim}{1}}_{g}❘{\hat{H}}_{dd}❘\overset{\sim}{1} \right\rangle\left. ❘{\overset{\sim}{1}}_{g} \right\rangle\left\langle {\overset{\sim}{1}}_{g}❘ \right.}}$ where $\left\langle {\overset{\sim}{0}}_{g}❘{\hat{H}}_{dd}❘{\overset{\sim}{0}}_{g} \right\rangle = {+ {\frac{\mu_{0}\mu_{N}{\gamma_{z}^{2}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}}{8\pi{r}^{3}\omega_{01}}\left\lbrack {{3\ln g_{vx}{\hat{I}}_{x}} + {3{mng}_{vx}{\hat{I}}_{y}} + {\left( {{3n^{2}} - 1} \right)g_{vz}{\hat{I}}_{z}}} \right\rbrack}}$ $\left\langle {\overset{\sim}{1}}_{g}❘{\hat{H}}_{dd}❘{\overset{\sim}{1}}_{g} \right\rangle = {- {{\frac{\mu_{0}\mu_{N}{\gamma_{z}^{2}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}}{8\pi{r}^{3}\omega_{01}}\left\lbrack {{3\ln g_{vx}{\hat{I}}_{x}} + {3{mng}_{vx}{\hat{I}}_{y}} + {\left( {{3n^{2}} - 1} \right)g_{vz}{\hat{I}}_{z}}} \right\rbrack}.}}$

Hence the effective interaction between the ¹⁷¹Yb qubit and the four register spins, Ĥ_(int)=Σ_(i∈register){tilde over (Ĥ)}_(dd) ^((i)), can be described by

${\hat{H}}_{int} = {{{\hat{\overset{\sim}{S}}}_{z}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}{\sum\limits_{i \in {register}}\left( {{J_{x}^{(i)}{\hat{I}}_{x}^{(i)}} + {J_{y}^{(i)}{\hat{I}}_{y}^{(i)}} + {J_{z}^{(i)}I_{z}^{(i)}}} \right)}}$ with $I_{x}^{(i)} = \frac{{- 3}\mu_{0}\mu_{N}\gamma_{z}^{2}g_{vx}l^{(i)}n^{(i)}}{4{\pi\left( r^{(i)} \right)}^{3}\omega_{01}}$ $I_{y}^{(i)} = \frac{{- 3}\mu_{0}\mu_{N}\gamma_{z}^{2}g_{vx}m^{(i)}n^{(i)}}{4{\pi\left( r^{(i)} \right)}^{3}\omega_{01}}$ $J_{z}^{(i)} = \frac{{- 3}\mu_{0}\mu_{N}\gamma_{z}^{2}{g_{vz}\left( {{3\left( n^{(i)} \right)^{2}} - 1} \right)}}{4{\pi\left( r^{(i)} \right)}^{3}\omega_{01}}$ and ${\overset{\hat{}}{\overset{\sim}{S}}}_{z} = {\frac{1}{2}{\left( {{\left. ❘{\overset{\sim}{1}}_{g} \right\rangle\left\langle {\overset{\sim}{1}}_{g}❘ \right.} - {\left. ❘{\overset{\sim}{0}}_{g} \right\rangle\left\langle {\overset{\sim}{0}}_{g}❘ \right.}} \right).}}$

Finally, local basis transformations of each ⁵¹V ion are performed to further simplify the Hamiltonian form. Specifically, we apply the following unitary rotation:

$\left. {\hat{H}}_{int}\rightarrow{U{\hat{H}}_{int}U^{\dagger}} \right.{{U = {\prod\limits_{j \in {register}}{\exp\left\lbrack {i\theta^{(j)}{\hat{I}}_{z}^{(j)}} \right\rbrack}}},{{{where}\theta^{(j)}} = {\tan^{- 1}\left( {m^{(j)}/I^{(j)}} \right)}},{{which}{leads}{to}}}{{\hat{H}}_{int} = {{{\hat{\overset{\sim}{S}}}_{z}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}{\sum\limits_{i \in {register}}\left\lbrack {{a_{x}{\hat{I}}_{x}^{(i)}} + {a_{z}{\hat{I}}_{z}^{(i)}}} \right\rbrack}}}{{{with}a_{x}} = {{\sqrt{\left( J_{x}^{(i)} \right)^{2} + \left( J_{y}^{(i)} \right)^{2}}{and}a_{z}} = {J_{z}^{(i)}.}}}$

Note that the coupling coefficients a_(x) and a_(z) are homogeneous (i.e. independent of site index i) since the four register spins are equidistant from the central ¹⁷¹Yb and have directional cosine factors with equal magnitude.

The same result can also be derived using the Schrieffer-Wolff transformation [4, 5], where the interaction Hamiltonian obtained here corresponds to the dominant second-order perturbation terms. Hereafter notation can be simplified by using |0_(g)

and |1_(g)

without tildes to represent the weakly perturbed eigenstates in the presence of any small magnetic field.

(vi) Full System Hamiltonian

Combining the various energy and interaction terms, the full system Hamiltonian (in a ¹⁷¹Yb frame rotating at ω₀₁/2π=675 MHz) becomes:

${\hat{H}}_{full} = {{\frac{{\gamma_{z}^{2}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right)}^{2}}{2\omega_{01}}{\overset{\hat{}}{\overset{\sim}{S}}}_{z}} + {\sum\limits_{i \in {register}}{{Q\left( {\hat{I}}_{z}^{(i)} \right)}^{2}{+ {\overset{\hat{}}{\overset{\sim}{S}}}_{z}}\left( {B_{z}^{OH} + {B^{RF}(t)}} \right){\sum\limits_{i \in {register}}{\left\lbrack {{a_{x}{\hat{I}}_{x}^{(i)}} + {a_{z}{\hat{I}}_{z}^{(i)}}} \right\rbrack.}}}}}$

b. Randomised Benchmarking and Dynamical Decoupling

High fidelity control of the ¹⁷¹Yb|0_(g)

↔|1_(g)

transition is essential for implementing the ZenPol sequence and enabling coherent ¹⁷¹Yb⁵¹V interactions. For example, a single swap operation realised by the ZenPol sequence contains 120 local ¹⁷¹Yb gates. Single qubit gate fidelity can be characterized using randomised benchmarking [6], which provides a value independent from state preparation or measurement (SPAM) errors. We apply randomly sampled single qubit Clifford gates constructed using π and π/2 rotations around the x and y directions followed by the single-gate inverse operation (FIG. 6 a ). When the number of gates, M_(gate), increases, the sequence error accumulates and the probability of returning to the initial |0_(g)

state reduces according to an exponential decay:

P=0.5+P _(n) d ^(M) ^(gate) .

When ensemble-averaged over a sufficiently large number of random gate sets (in our case 100), f=½(1+d) becomes a reliable estimate of the average single-qubit gate fidelity. Measurement results are presented in FIG. 6 a , leading to an extracted average single qubit gate fidelity of f=0.99975±0.00004.

The T₂ coherence time of the qubit transition is measured using an XY-8 dynamical decoupling sequence [7]. Specifically, we work with a fixed inter-pulse separation of 5.6 μs and measure the coherence time by varying the number of decoupling periods, M^(r) (FIG. 6 b ). An exponential decay 1/e with time constant T₂=16±2 ms is measured. This measurement uses the same method as in [8], however, we observe a factor of three improvement in coherence due to the improved microwave setup leading to correspondingly increased π gate fidelities.

c. Extra Register Detail

In this section additional technical details are provided related to the single excitation states used to store quantum information on the ⁵¹V spins.

The general form for the engineered spin-exchange interaction is:

${{\hat{H}}_{avg} = {B^{RF}{\sum\limits_{i \in {register}}{b_{({k,\omega_{j}})}^{(i)}\left( {{{\overset{\hat{}}{\overset{\sim}{S}}}_{+}{\overset{\hat{}}{\overset{\sim}{I}}}_{-}^{(i)}} + {{\overset{\hat{}}{\overset{\sim}{S}}}_{-}{\overset{\hat{}}{\overset{\sim}{I}}}_{+}^{(i)}}} \right)}}}},$

where, B^(RF) is the square-wave RF magnetic field amplitude, b_((k,ω) _(j) ₎ ^((i)) is a k-dependent prefactor for transitions ω_(i) of the i^(th) register spin, {tilde over (Î)}₊=|↑

↓|, {tilde over (Î)}⁻=|↓

(↑|, are raising and lowering operators in an effective nuclear spin-1/2 manifold and {tilde over (Ŝ)}₊=|1_(g)

0_(g)|, {tilde over (Ŝ)}⁻=|0_(g)

1_(g)| are raising and lowering operators for the ¹⁷¹Yb qubit. Note, in this section, we do not assume homogeneous coupling to the register spins, hence the b_((kωj)) ^((i)) coefficients depend on the register site index i. In addition, we consider an arbitrary number of register spins, N, that are spectrally indistinguishable.

When the ¹⁷¹Yb is initialised in |1_(g)

and the ⁵¹V register spins are polarised in |0_(v)=|↓

^(⊗N), this interaction leads to the following spin-exchange evolution [9]:

|ψ(t)

=|1_(g)

|0_(v)

cos(J _(ex) t/2)−i|0_(g)

|1_(v)

sin(J _(ex) t/2)

where the spin-exchange frequency is given by:

$I_{ex} = {2B^{rf}\sqrt{\sum\limits_{i \in {register}}{❘b_{({k,\omega_{j}})}^{(i)}❘}^{2}}}$

and the resulting single-spin excited state generated by this interaction is:

$\left. ❘1_{v} \right\rangle = {\frac{1}{\Sigma_{i}{❘b_{({k,{\omega}_{j}})}^{(i)}❘}^{2}}{\sum\limits_{i \in {register}}\left. b_{({k,{\omega}_{j}})}^{(i)}❘\downarrow\ldots\uparrow{}_{(i)}\ldots\downarrow{\rangle.} \right.}}$

Based on the results presented in FIG. 2 d and Supplementary Information Section VIII we postulate that for our system the register consists of the second nearest shell of four homogeneously coupled ⁵¹V ions. In this case we recover the expressions presented in the main text, namely, the single-spin excitation in the register realises an entangled four-body W-state, |W_(v)

, as depicted in FIG. 1 c :

$\left. ❘1_{v} \right\rangle = {\left. ❘W_{v} \right\rangle = \frac{\text{↑↓↓↓〉+↓↑↓↓〉+↓↓↑↓〉+↓↓↓↑〉}}{2}}$

and the spin-exchange rate is given by J_(ex)=4B^(RF)b_((k,ω) _(j) ₎. In general, for N homogeneously coupled register spins, we expect that the spin-exchange rate is enhanced by a factor of √{square root over (N)}, leading to faster swap gate operation.

In this protocol it is possible to transfer a second spin excitation to the register. More specifically, the spin-preserving exchange interaction, {tilde over (Ŝ)}−{tilde over (Î)}₊+{tilde over (Ŝ)}+{tilde over (Î)}⁻, couples the state |1_(g)

|W_(v)

to |0_(g)

|2_(v)

, where |2_(v)

is a ⁵¹V state with two spins in |↑

. To avoid undesired excitation to states outside of the effective {0_(v)

, |W_(v)

} manifold, the ¹⁷¹Yb qubit in |0_(g)

is always prepared before retrieving stored states from the ⁵¹V register. Hence the swap gate realised by this interaction operates on a limited basis of states.

Utilising the dense, lattice nuclear spins ensures near identical registers for all ¹⁷¹Yb ions. FIG. 14 shows ZenPol spectra near the ω_(c) transition, collectively enhanced spin-exchange oscillations and motionally-narrowed T₂* times for three ⁵¹V registers coupled to three different ¹⁷¹Yb ions. The ¹⁷¹Yb optical and microwave frequencies were re-calibrated for each ion, however, all aspects of the experimental sequences related to register control and readout were identical.

d. Simulation

We simulate our coupled spin system using the effective Hamiltonian derived in Supplementary Information, however we add three additional terms:

1 Nuclear Zeeman interactions of the ⁵¹V register spins with the Overhauser field from the bath: Since the energy levels are quantised along the z-axis, magnetic fluctuations along the z-direction dominate, which can be captured by the following Hamiltonian

$\text{?} = {\sum\limits_{i \in {register}}{\mu_{N}g_{vz}{B_{z}^{OH}\left( r_{i} \right)}{\hat{I}}_{z}^{(i)}}}$ ?indicates text missing or illegible when filed

where B_(π) ^(OH)(r_(i)) is the z-component of the Overhauser field evaluated at the position of the i^(th) register ion, r_(i).

2 Nuclear magnetic dipole-dipole interactions of the register spins:

$\text{?} = {\underset{i < j}{\sum\limits_{i,{j \in {register}}}}{\frac{\mu_{0}}{4\pi}\left\lbrack {\frac{\mu_{v}^{(i)} \cdot \mu_{v}^{(j)}}{r_{ij}^{3}} - \frac{3\left( {\mu_{v}^{(i)} \cdot r_{ij}} \right)\left( {\mu_{v}^{(j)} \cdot r_{ij}} \right)}{r_{ij}^{5}}} \right\rbrack}}$ ?indicates text missing or illegible when filed

with r_(ij) the displacement vector between ⁵¹V register spins at site i and j. 3. ¹⁷¹Yb-enhanced register spin-spin interactions: These terms are derived by considering second-order perturbations using the Schrieffer-Wolff transformation [4,5]. For example, the dominant Ising-type terms take the form

${{\hat{H}}_{edd} = {\sum\limits_{i,{j \in {register}}}{{\frac{1}{2\omega_{01}}\left\lbrack {\left( {{3n^{2}} - 1} \right)\frac{\mu_{0}\mu_{N}\gamma_{z}g_{vz}}{4\pi r^{3}}} \right\rbrack}^{2}{\hat{S}}_{z}{\hat{I}}_{z}^{(i)}{\hat{I}}_{z}^{(j)}}}},$

where r and n are the magnitude and c-direction cosine of the ¹⁷¹Yb ⁵¹V register ion displacement vector. However, we note that the ZenPol sequence cancels these interactions to first order.

By simulating ¹⁷¹Yb Ramsey coherence times, g_(vz)≈1.6 is extracted. Estimation of the bare ⁵¹V coherence time indicates a potential discrepancy in this value g_(vz) by up to 25%, discussed further in Supplementary Information, however, this has a negligible impact on the ZenPol sequence simulations. An estimate for g_(vx)≈0.6 is obtained by calibrating the RF field amplitude and comparing with the experimental results of direct ⁵¹V spin driving in FIG. 11 .

The nuclear Overhauser field B_(π) ^(OH) is computed according to equation (S8) by randomly sampling the bath states for each Monte-Carlo simulation repetition. A simple model of the bath dynamics is included by incorporating stochastic jumps of the bath spins on magnetic-dipole allowed transitions.

The register spin dynamics are simulated in a reduced Hilbert space by considering only the ω_(c) manifold. This enables fast simulation of all four register spins plus the ¹⁷¹Yb qubit transition (Hilbert space with dimension 32). Imperfect polarisation of the ⁵¹V register into |↓

=|±7/2

is categorised into two distinct types:

1 Imperfect polarisation within the ω_(c) transition i.e. a small residual population ϵ₁ in |↑

=|5/2

.

2 Imperfect polarisation outside the ω_(c) manifold i.e. a small residual population ϵ₂ in |±1/2

and |±3/2

.

This leads to a |↓

population of 1−ϵ₁−ϵ₂. Incomplete polarization is incoporated by sampling different register initial states for each Monte-Carlo repetition. For case 1, this involves occasionally initialising a given ⁵¹V ion into |↑

, while for case 2 this involves reducing the Hilbert space dimension by removing the ⁵¹V ion from the simulation. Finite pulse duration effects are taken into account by modeling the ZenPol sequence using 25 nsπ/2 and 50 nsπ pulses (FIG. 10 a ).

As shown in Extended Data FIG. 6 d , the spin-exchange oscillations from numerical simulation (red dashed line) exhibit slower decay than the measured experimental results (red markers). A phenomenological exponential decay envelope, ce^(−M/τ) ^(M) , is added to the simulation results where c and τ_(M) are free parameters, and M is the ZenPol sequence period. The additional decay could be caused by heating due to the RF field, excess ¹⁷¹Yb dephasing or additional register spin interactions which we haven't considered here. This model is fitted by optimising multiple parameters: ϵ₁,ϵ₂,B^(RF),c and τ_(M). The resulting values of ϵ₁ and ϵ₂ are 0.12 and 0.04, respectively, indicating ≈84% polarisation into |↓

; the RF magnetic field amplitude is B^(RF)≈1.6G and the phenomenological exponential decay parameters are c=0.8 and τ_(M)=90 leading to a close fit with the experimental results (red solid line, FIG. 2 c and FIG. 10 d ). Additional simulation results following this methodology with varying B_(RF) and τ are presented in FIG. 10 .

Finally, the results are modeled with a single-spin excitation in the ω_(c)-manifold by including the |±3/2

level in the simulation (FIG. 2 d and Supplementary Information). The initial state used in this simulation is partially polarised between the |±3/2

level with population 1−ϵ and the |±1/2

level with population ϵ. We use the same value of B^(RF)=1.6G as in FIG. 2 c , and optimise the polarisation level leading to 1−ϵ=0.8. The close correspondence between the measured and simulated oscillation profiles suggests that the register consists of the second shell of four homogeneously coupled ⁵¹V ions.

e. Hartmann Hahn Spectroscopy

In addition to the ZenPol spectra discussed above, Hartmann-Hahn (HH) double resonance [10] is used to perform spectroscopy of the nuclear spin environment. This method enables spin exchange between two systems with different transition frequencies by resonantly driving a qubit with a Rabi frequency that matches the energy level splitting of the environmental nuclear spins. In our case, we resonantly drive the ¹⁷¹Yb at 675 MHz to generate a pair of dressed states

$\left. ❘ \pm \right\rangle = {\frac{1}{\sqrt{2}}\left( {\left. ❘0_{g} \right\rangle \pm \left. i❘1_{g} \right\rangle} \right)}$

with splitting Ω which we sweep over a range ≈2π×(0−2.3) MHz (FIG. 7 ). The ¹⁷¹Yb qubit is initialised into the |−

dressed state by a π/2 pulse preceding the driving period. If resonant with a nuclear spin transition, the ¹⁷¹Yb qubit undergoes spin exchange at a rate dictated by the interaction strength. Finally we read out the ¹⁷¹Yb|+

dressed state population to determine whether spin exchange has occurred. FIG. 7 b shows experimental results of Hartmann-Hahn spectroscopy where we vary both the HH drive Rabi frequency (Ω) and also the HH pulse duration (t). The counts plotted on the colour-bar are proportional to the |+

dressed state population. Three clear resonances are found at evenly spaced pulse amplitudes 0.15, 0.30 and 0.45 corresponding to the ω_(a), ω_(b) and ω_(c) ⁵¹V transitions; notably, unlike ZenPol, the HH sequence only has one harmonic leading to a single resonant interaction per transition. Also note the lack of oscillations when varying the pulse duration, t, on resonance with either of the three transitions: this is because the spin exchange is driven by the randomised, Overhauser field induced ¹⁷¹Yb dipole moment. For this reason, the HH sequence cannot be used to generate the coherent exchange interaction necessary to realise a swap gate for our system. In the case of no driving (Ω=0), the signal rapidly saturates as t increases as a result of Ramsey dephasing of the initial state. However, as Ω exceeds the ¹⁷¹Yb spin linewidth (˜50 kHz[8]), this effect diminishes due to the emergence of spin-locking effects and consequently leads to an increased saturation timescale when not resonant with the ⁵¹V transitions. The resolution of this measurement is also limited by the ¹⁷¹Yb spin linewidth, and we therefore cannot resolve the split-resonance structure observed in the ZenPol spectra. The results agree well with simulations (FIG. 7 c ) verifying that interactions with the ⁵¹V quadrupolar structure dominate these measurements.

f. Polarisation of Multi-Level Register Nuclear Spins

Polarisation dynamics are explored using the PROPI method (polarisation readout by polarisation inversion) [11]. This sequence uses the back-action of the ⁵¹V spins on the ¹⁷¹Yb to measure the register polarisation after successive ZenPol polarisation cycles. For instance, when polarising into |↑

=|±5/2

on the ω_(c) transition, the ¹⁷¹Yb is initialised into |1_(g)

and undergoes spin exchange with any ⁵¹V population in |↓

=|±7/2

. The ¹⁷¹Yb|0_(g)

population after interaction is therefore related to the residual ⁵¹V|↓

population. As presented in FIG. 9 a , the ¹⁷¹Yb population is measured after each of 20 consecutive polarisation cycles and a saturation is observed after 10 cycles, indicating that the ¹⁷¹Yb polarisation has been transferred to the ⁵¹V register. The high-contrast signal obtained in this measurement is enabled by alternating the ⁵¹V polarisation direction, i.e. periods of polarisation into |↑

are interleaved with periods of polarisation into |↓

. This mitigates the need to wait for slow register thermalisation (T₁ ⁽⁰⁾≈0.5 s, see Supplementary Information Section X) between consecutive experiment repetitions. These measurements are repeated with ZenPol sequences on the ω_(b) transition, demonstrating similar levels of polarisation saturation after approximately 10 cycles (FIG. 9 b ).

We also demonstrate the effect of incomplete register polarization on the spin-exchange oscillation by varying the number of polarisation cycles on the ω_(b) and ω_(c) transitions before each experiment (FIG. 9 c ). As expected, the coherent spin-exchange oscillations emerge as an increasing number of polarisation cycles are applied.

These results inform the design of polarisation sequences used in subsequent single-spin excitation experiments where 40 polarisation cycles interleaved between the ω_(b) and ω_(c) transitions are sufficient to polarise the register into |0_(v)

=|↓↓↓↓

. Based on simulations discussed in Supplementary Information, we estimate this protocol achieves ≈84% polarisation into the |0_(v)

state. Note the ZenPol sequence is not used to directly polarise the ω_(a) transition due to spectral overlap with ω_(b) and ω_(c) (FIG. 2 b ). We postulate that the high degree of polarisation can still be achieved even in the absence of direct ω_(a) transition control due to two factors:

1 The thermalisation timescale of the ω_(a) transition is significantly shorter than the interrogation time. Specifically, our experiments typically run for several minutes whereas the ω_(a) thermalisation rate is likely similar to T₁ ⁽⁰⁾=0.54 s. Thus, undesired population in the |±1/2

level can still pumped to |±7/2

once it relaxes to |±3/2

.

2 Once successfully initialised into the ω_(c) manifold the probability of shelving into the |±1/2

level is small as it necessitates two consecutive decays on the ω_(b) and ω_(a) transitions, both of which are considerably slower than our experiment/polarisation repetition rate (20 ms).

We tried to improve the polarisation fidelity by incorporating direct driving on the ω_(a) transition during the polarisation protocol. This leads to fast population exchange between ±1/2

and |±3/2

, however, there was no improvement to the contrast of the resulting spin exchange oscillations thereby indicating that shelving into |±1/2

is not a limiting factor in our experiments.

g. Analysis of Spin Exchange Dynamics

In this section, an analysis of the spin exchange dynamics on the ω_(c) register transition is presented. The spin-exchange measurements in FIG. 2 c are measured at a fixed ZenPol period of 2τ=5.048 μs leading to resonant interactions with the 991 kHz ω_(c)-transition. However, analogous to the Rabi oscillations in a two-level system, the oscillation frequency and contrast of these spin transfer oscillations also depend on the detuning of the ZenPol sequence relative to the ⁵¹V transition. Specifically, we expect the following relations:

${{J_{ex}(\delta)} = \sqrt{{J_{ex}(0)}^{2} + \delta^{2}}}{{C(\delta)} = \frac{{J_{ex}(0)}^{2}}{{J_{ex}(0)}^{2} + \delta^{2}}}$

Here J_(ex) and C are the spin-exchange frequency and oscillation contrast, respectively, and δ is the detuning of the ZenPol sequence resonance relative to a target nuclear spin transition. The register is polarized into |0_(v)

and FIG. 10C shows measurement of the frequency detuning dependence of the spin-exchange oscillations. These results agree well with the corresponding simulations shown in FIG. 10 c.

Control of the spin exchange frequency by varying the RF magnetic field amplitude (B^(RF)) is also demonstrated. FIG. 10D shows the spin-exchange dynamics for four different values of B^(RF)=0.8G, 1.2G, 1.6G and 2.0G. The inset in FIG. 2C plots extracted spin exchange frequencies J_(ex) for a range of different B^(RF) demonstrating linear dependence as expected and leading to accurate control of the engineered interaction strength (see First Example for details).

h. Single Excitation in ω_(c) Manifold

The ability to shelve populations in different quadrupole levels enables the operation of the ⁵¹V register with an alternative set of many-body states: |0_(v)′

and |1_(v)′

. For this experiment the ⁵¹V spins are polarized down the energy ladder on the ω_(b) and ω_(c) transitions leading to polarisation primarily into the |±3/2

level, with a small residual population in |±1/2

. For the purpose of this analysis, perfect polarisation into |±3/2

is assumed, however ω_(a) transition polarisation/addressability would be required for this.

The register |1_(v)′

state is prepared by injecting a single spin excitation on the ω_(b) transition (i.e. from |±3/2

→|↑

=|±5/2

), this is achieved using the corresponding ZenPol resonance at ω_(b),k=3:

$\left. \left. {\left. {\left. {\left. {\left. {❘1_{v}^{\prime}} \right\rangle = {\frac{1}{2}\left( {❘\left. \uparrow,\frac{3}{2},\frac{3}{2},\frac{3}{2} \right.} \right.}} \right\rangle + {❘\left. \frac{3}{2},\uparrow,\frac{3}{2},\frac{3}{2} \right.}} \right\rangle + {❘\left. \frac{3}{2},\frac{3}{2},\uparrow,\frac{3}{2} \right.}} \right\rangle + {❘\left. \frac{3}{2},\frac{3}{2},\frac{3}{2},\uparrow \right.}} \right\rangle \right)$

Here the ± sign is omitted in the state label for simplicity. Subsequently, the ¹⁷¹Yb in |0_(g)

and induce a spin exchange oscillation between |↑

and |↓

=|±7/2

via a ZenPol sequence resonant with the ω_(c) transition. The resulting time evolution is given by

$\left. ❘{\psi(t)} \right\rangle = {{\left. ❘0_{g} \right\rangle\left. ❘1_{v}^{\prime} \right\rangle\cos\left( \frac{J_{ex}^{\prime}t}{2} \right)} - {\left. i❘1_{g} \right\rangle\left. ❘0_{v}^{\prime} \right\rangle\sin\left( \frac{J_{ex}^{\prime}t}{2} \right)}}$ where $\left. \left. {\left. {\left. {\left. {\left. {❘0_{v}^{\prime}} \right\rangle = {\frac{1}{2}\left( {❘\left. \downarrow,\frac{3}{2},\frac{3}{2},\frac{3}{2} \right.} \right.}} \right\rangle + {❘\left. \frac{3}{2},\downarrow,\frac{3}{2},\frac{3}{2} \right.}} \right\rangle + {❘\left. \frac{3}{2},\frac{3}{2},\downarrow,\frac{3}{2} \right.}} \right\rangle + {❘\left. \frac{3}{2},\frac{3}{2},\frac{3}{2},\downarrow \right.}} \right\rangle \right)$

and J_(ex)′=2b_((k,ω) _(c) ₎B^(RF) with k=5. Notice that the spin-exchange oscillation rate, J_(ex)′, no longer has a √{square root over (N)} rate enhancement, this is because every ket in the |1_(v)′

and |0_(v)′

states contains only a single spin in the ω_(c)-transition manifold. Using this manifold for information storage would have several benefits. For instance, direct microwave driving of the register ω_(c)-transition would lead to Rabi oscillation between |0_(v)′

and |1_(v)′

and could therefore be used to realise local gates in this basis. Additionally, a second spin excitation is not allowed in this scheme, therefore the ZenPol sequence reproduces a complete two-qubit swap gate regardless of the ¹⁷¹Yb state. For these reasons, we believe that there may be some advantages to working with the {|0_(v)′

, |1_(v)′

} manifold if the state initialisation fidelity into |±3/2

can be improved via direct ω_(a) transition polarisation.

i. T₂* Coherence Discussion

Here we provide detailed discussions regarding the ⁵¹V register coherence decay processes described in the main text. There are two magnetic interactions which limit the T₂* dephasing timescale: (1) the direct nuclear Zeeman interaction of each register spin with the Overhauser field (equation (S20)) and (2) a contribution from the ¹⁷¹Yb Knight field [12]. In the latter case, the bath-induced ¹⁷¹Yb dipole moment generates a randomly fluctuating magnetic field at each ⁵¹V ion, the Knight field, which is described by

? = ±g_(vz)μ_(N)B_(z)^(OH)A_(z)Î_(z) with $A_{z} = \frac{\mu_{8}{\gamma_{z}^{2}\left( {1 - {3n^{2}}} \right)}}{8\pi r^{3}\text{?}}$ ?indicates text missing or illegible when filed

Here, the + and − cases in equation (S28) correspond to ¹⁷¹Yb in |1_(g)

and |0_(g)

, respectively. The constants are defined in Supplementary Information Section. We note that A_(z) corresponds to an effective local field amplification factor with value A_(z)≈3.1 for the register spins. We define the ¹⁷¹Yb Knight field to be ±A_(π)B_(π) ^(OH).

By applying periodic π pulses to the ¹⁷¹Yb, we flip its state between |0_(g)

and |1_(g)

, thereby switching the sign of the Knight field. This leads to the cancellation of V phase accumulation between successive free evolution periods, resulting in a longer coherence time. We numerically simulate the register coherence times using the method outlined in Supplementary Information Section IV. When limited by the ¹⁷¹Yb Knight field, simulation yields a Gaussian decay with a 1/e coherence time of 33 μs (equivalent to experimental results in FIG. 3 a ). We also predict an upper bound for the coherence time when decoupled from the ¹⁷¹Yb Knight field by turning off Hamiltonian terms associated with equation (S28), yielding an extended Gaussian decay of 417 μs (equivalent to experimental results in FIG. 3 b ). These simulated values are consistent with the corresponding experimental results (58±4 μs and 225±9 μs respectively) to within a factor of two. We note that this could indicate an error in our estimation of g_(vz) by up to 25%, potentially caused by a small discrepancy in the position of the two ⁵¹V bath spins closest to ¹⁷¹Yb. Further analysis of these parameters is left for future work.

j T₁ Lifetime Discussion

We measure the population decay of both the |0_(v)

and |W_(v)

states (timescales T₁ ⁽⁰⁾ and T₁ ^((W)) respectively) by preparing the ⁵¹V register in the appropriate state and waiting for a variable time, t, before swapping to the ¹⁷¹Yb for readout.

The |0_(v)

state exhibits slow exponential decay with 1/e time constant T₁ ⁽⁰⁾=0.54±0.08 s (FIG. 12 b ). There are two contributions which could be limiting this decay:

1 Resonant population exchange between the register spins and unpolarised frozen-core ‘dark spins’. For instance, the two nearest ⁵¹V ions (ions 1 and 2 in the table in Supplementary Information Section) may interact resonantly with the neighbouring register spins. However, we cannot detect or polarise these dark spins since they only interact with the ¹⁷¹Yb via Ising-like Ŝ_(π)Î_(π) terms.

2 Off-resonant population exchange between the register and detuned unpolarised bath spins.

As for the |W_(v)

state, it exhibits a Gaussian decay with a much faster 1/e time constant of T₁ ^((W))=39.5±1.3 μs (FIG. 12 a ). This can be explained by considering the effect of dephasing on the register spins. Specifically, the |W_(v)

state which our ¹⁷¹Yb qubit interacts with is given as

|W _(v)

=½(↑↓↓↓

+|↓↑↓↓

+|↓↓|↓

|↓↓↓↑

).

Crucially, there are three additional orthogonal states required to span the ⁵¹V register single excitation subspace:

|α_(v)

=½(|↑↓↓↓

+|↓↑↓↓

−|↓↓↑↓

−|↓↓↓↑

)

|β_(v)

=½(|↑↓↓↓

−|↓↑↓↓

+|↓↓↑↓

−|↓↓↓↑

)

|γ_(v)

=½(|↑↓↓

−|↓↑↓↓

−|↓↓↑↓

+|↓↓↓↑

)

We assume uncorrelated noise at each of the four ⁵¹V spins and apply a pure-dephasing master equation model. In the single excitation subspace, this becomes:

? = 2Γ[𝒟(↑↓↓↓)(↑↓↓)+𝒟(↓↑↓↓)(↓↑↓↓)+𝒟(↓↓↑↓)(↓↓↑↓)+𝒟(↓↓↓↑)(↓↓↑)]ρ ?indicates text missing or illegible when filed

where the dephasing channel (Lindbladian) is given by)

(â)ρ=âρâ ^(†)−½{â ^(†) â,ρ}

and Γ is the dephasing rate on the ω_(c) transition of a single ⁵¹V spin. We solve this equation for different initial states ρ(0). When ρ(0)=|0_(v)

0_(v)|, dephasing does not contribute to T₁ ⁽⁰⁾, i.e. ρ(t)=ρ(0). However, when ρ(0)=|W_(v)

W_(v)| the state evolves according to

ρ(t)=|W _(v)

W _(v) |e ^(−2Γt)+¼(1−e ^(−2Γt))1^((SEM))

where 1^((SEM)) is the single excitation manifold identity operator:

1^((SEM)) =|W _(v)

W _(v)|+|α_(v)

α_(v)|+|β_(v)

β_(v)|+|γ_(v)

γ_(v)|

i.e. dephasing leads to decay of |W_(v)

into 1^((SEM)) at rate 2Γ. For completeness we also consider the decay of the off-diagonal coherence term ρ₀₁=

0_(v)|ρ|W_(v)

and find that

ρ₀₁(t)=ρ₀₁(0)e ^(−Γt).

Essentially, the pure dephasing model predicts T₂*=2T₁ ^((W)) for our system.

We verify that dephasing is the main source of |W_(v)

population decay by demonstrating lifetime extension using the same motional narrowing approach employed to improve the coherence time. Specifically, during the wait time, we apply a series of π pulses to the ¹⁷¹Yb separated by 6 μs leading to an extended lifetime of T₁ ^((W))=127±8 μs (FIG. 12 a ). We note that both the bare and motionally-narrowed T₁ ^((W)) and T₂* times are close to the T₂*=2T₁ ^((W)) limit identified above. We further extend the T₁ ^((W)) lifetime to 640±20 μs using two ⁵¹Vπ pulses applied during the wait time, thereby achieving dynamical decoupling from the nuclear Overhauser field (equivalent to the results in FIG. 3 c ).

Finally we note that if T₁ ^((W)) is limited by the ¹⁷¹Yb Knight field as a common noise source, there may be some discrepancy in the predictions of this model due to a high degree of noise correlation between the four ⁵¹V register spins arising from lattice symmetry. However, when performing motional narrowing we decouple the ¹⁷¹Yb Knight field and are likely limited by the, considerably less correlated, local Overhauser field.

k. Parity Oscillations and Coherence

Here we derive an expression for the ¹⁷¹Yb ⁵¹V Bell-state coherence ρ₀₁=

1_(g)0_(v)|ρ|0_(g)W_(v)

in terms of the parity oscillation contrast with a correction factor. In particular, when reading out this coherence, we apply a √{square root over (swap)} gate which maps |Ψ⁺

=½(|1_(g)0_(v)

−i|0_(g)W_(v)

) to |0_(g)W_(v)

and |Ψ⁻

=½(|1_(g)0_(v)

+i|0_(g)W_(v)

) to |1_(g)0_(v)

. Note that reading out the ¹⁷¹Yb state is sufficient to distinguish the |Ψ⁺

and |Ψ⁻

states in this measurement. We can account for the readout fidelity of the |Ψ^(±)

states by using a √{square root over (

_(,1))} factor (Methods), i.e. if the state |Ψ⁺

(|Ψ⁻

) is perfectly prepared, ¹⁷¹Yb will be measured in state |0_(g)

(|1_(g)

) with probability ½(1+√{square root over (

_(,1))}). To span the ¹⁷¹Yb-⁵¹V Hilbert space, we also need to consider the effect of the readout √{square root over (swap)} gate when the system is initialised into the other two states: |1_(g)W_(v)

or |0_(g)0_(v)

. To this end, we assign imperfect readout probabilities of q₁₁ and q₀₀ for |1_(g)W_(v)

and |0_(g)0_(v)

, respectively. Specifically, we can represent the dependence of the parity readout on the input state using the following matrix relation:

$\begin{pmatrix} p_{1,{Yb}} \\ p_{0,{Yb}} \end{pmatrix} = {\mathcal{M}_{swap}{\mathcal{M}_{wait}\begin{pmatrix} p_{11} \\ p_{\Psi^{+}} \\ p_{\Psi^{-}} \\ p_{00} \end{pmatrix}}}$ with $\mathcal{M}_{swap} = \begin{pmatrix} q_{11} & {\frac{1}{2}\left( {1 - \sqrt{\mathcal{F}_{{sw},1}}} \right)} & {\frac{1}{2}\left( {1 + \sqrt{\mathcal{F}_{{sw},1}}} \right)} & {1 - q_{00}} \\ {1 - q_{11}} & {\frac{1}{2}\left( {1 + \sqrt{\mathcal{F}_{{sw},1}}} \right)} & {\frac{1}{2}\left( {1 - \sqrt{\mathcal{F}_{{sw},1}}} \right)} & q_{00} \end{pmatrix}$ $\mathcal{M}_{wait} = {\begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & {\cos^{2}\left( {\text{?}t/2} \right)} & {\sin^{2}\left( {\text{?}t/2} \right)} & 0 \\ 0 & {\sin^{2}\left( {\text{?}t/2} \right)} & {\cos^{2}\left( {\text{?}t/2} \right)} & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix}.}$ ?indicates text missing or illegible when filed

Here p_(1,Yb) and p_(0,Yb) are the probabilities of measuring the ¹⁷¹Yb qubit in |1_(g)

and |0_(g)

, respectively, and p_(Ψ) _(±) =

Ψ^(±)|ρ|Ψ^(±)

are the probabilities of being in the |Ψ^(±)

Bell states. The contrast C_(parity) of the parity oscillation between |105 ⁺

and |Ψ⁻

is extracted by measuring the difference in the ¹⁷¹Yb|1_(g)

populations measured at t=0 and t=π/ω_(c), allowing us to estimate the Bell state coherence as |ρ₀₁|=C_(parity)/2√{square root over (

_(,1))}. This implies that uncorrected and corrected Bell state coherence values differ by a factor of √{square root over (

_(,1))}=0.72. Using the results presented in FIG. 4 b we obtain corrected and uncorrected estimates for |ρ₀₁| of 0.35.2±0.004 and 0.254±0.003 respectively.

1. Bell State Fidelity Error Analysis

To extract the Bell state fidelity and uncertainty, we perform a maximum likelihood analysis of the population and parity oscillation measurements, adopting a similar approach as in [13]. The population measurement involves a series of n experiments with outcomes distributed between the four population states: {n₀₀, n₀₁, n₁₀, n₁₁} where n=₀₀+n₀₁+n₁₀+n₁₁. The likelihood function for the uncorrected populations, {p₀₀, p₀₁, p₁₀, p₁₁} has multinomial form:

${\mathcal{L}\left( \left\{ p_{ij} \right\} \middle| \left\{ n_{ij} \right\} \right)} = {\frac{n!}{{n_{00}!}{n_{01}!}{n_{10}!}{n_{11}!}}p_{00}^{n_{00}}p_{01}^{n_{01}}p_{10}^{n_{10}}p_{11}^{n_{11}}}$

where we have assumed a prior uniform over the physical values of {p_(ij)}, i.e. 0≤p_(ij)≤1 and Σp_(ij)=1. The likelihood function for the corrected populations, {c₀₀, c₀₁, c₁₀, c₁₁}, is obtained by substituting equation (11) into equation (S35) and assuming a prior uniform over the physical values of {c_(ij)}, i.e. 0≤c_(ij)≤1 and Σc_(ij)=1. Corrected populations are obtained by maximising this likelihood function. The error for a specific population (say c₀₀) is obtained by marginalising

({c_(ij)}|{n_(ij)}) over the other three (c₀₁, c₁₀, c₁₁) and taking a 68% symmetric confidence interval.

We extract a likelihood function for the coherence by considering the following model:

y _(i)=0.5+√{square root over (

_(,i))}ρ₀₁ cos(ω_(c) t _(i))+ϵ_(i)

where {t_(i),y_(i)} are the parity oscillation data at the i^(th) point, ρ₀₁ is the corrected coherence,

_(sw,1) is the parity oscillation correction factor associated with the swap gate infidelity, and ϵ_(i) is the experimental error assumed to be normally distributed with μ=0 and unknown σ. The likelihood function is given by

${\mathcal{L}\left( {\rho_{01},\left. \sigma \middle| \left\{ {t_{i},y_{i}} \right\} \right.} \right)} = {\prod\limits_{i}{\frac{1}{\sqrt{2\pi}\sigma}{{\exp\left\lbrack \frac{\left( {y_{i} - 0.5 - {\sqrt{\mathcal{F}_{{sw},1}}\rho_{01}\cos\left( {\omega_{c}t_{i}} \right)}} \right)^{2}}{2\sigma^{2}} \right\rbrack}.}}}$

We obtain a likelihood for the corrected coherence,

(ρ₀₁|{t_(i),y_(i)}) by marginalising over σ.

The likelihood function for the fidelity is obtained by taking a product of the likelihood function for the populations with the likelihood function for the coherence and evaluating a contour integral at constant

, given as

${\mathcal{L}(\mathcal{F})} = {\int_{\mathcal{F}}{{\mathcal{L}\left( \left\{ c_{ij} \right\} \middle| \left\{ n_{ij} \right\} \right)}{\mathcal{L}\left( \rho_{01} \middle| \left\{ {t_{i},y_{i}} \right\} \right)}d\rho_{01}{\prod\limits_{ij}{{dc}_{ij}.}}}}$

The Bell state fidelity is extracted by maximising this likelihood and the error is evaluated as a symmetric 68% confidence interval.

Hardware Environment

FIG. 15 is an exemplary hardware and software environment 1500 (referred to as a computer-implemented system and/or computer-implemented method) used to implement one or more embodiments of the invention. The hardware and software environment includes a computer 1502 and may include peripherals. Computer 1502 may be a user/client computer, server computer, or may be a database computer. The computer 1502 comprises a hardware processor 1504A and/or a special purpose hardware processor 1504B (hereinafter alternatively collectively referred to as processor 1504) and a memory 1506, such as random access memory (RAM). The computer 1502 may be coupled to, and/or integrated with, other devices, including input/output (I/O) devices such as a keyboard 1514, a cursor control device 1516 (e.g., a mouse) a pointing device, pen and tablet, touch screen, multi-touch device, etc.) and a printer 1528. In one or more embodiments, computer 1502 may be coupled to, or may comprise, a portable or media viewing/listening device 1532. In yet another embodiment, the computer 1502 may comprise a multi-touch device, mobile phone, or other internet enabled device executing on various platforms and operating systems.

In one embodiment, the computer 1502 operates by the hardware processor 1504A performing instructions defined by the computer program 1510 under control of an operating system 1508. The computer program 1510 and/or the operating system 1508 may be stored in the memory 1506 and may interface with the user and/or other devices to accept input and commands and, based on such input and commands and the instructions defined by the computer program 1510 and operating system 1508, to provide output and results.

Output/results may be presented on the display 1522 or provided to another device for presentation or further processing or action. The image may be provided through a graphical user interface (GUI) module 1518. Although the GUI module 1518 is depicted as a separate module, the instructions performing the GUI functions can be resident or distributed in the operating system 1508, the computer program 1510, or implemented with special purpose memory and processors.

Some or all of the operations performed by the computer 1502 according to the computer program 1510 instructions may be implemented in a special purpose processor 1504B. In this embodiment, some or all of the computer program 1510 instructions may be implemented via firmware instructions stored in a read only memory (ROM), a programmable read only memory (PROM) or flash memory within the special purpose processor 1504B or in memory 1506. The special purpose processor 1504B may also be hardwired through circuit design to perform some or all of the operations to implement the present invention. Further, the special purpose processor 1504B may be a hybrid processor, which includes dedicated circuitry for performing a subset of functions, and other circuits for performing more general functions such as responding to computer program 1510 instructions. In one embodiment, the special purpose processor 1504B is an application specific integrated circuit (ASIC) or field programmable gate array (FPGA). In other examples, special purpose processor may comprise a graphics processing unit (GPU).

The computer 1502 may also implement a compiler 1512 that allows an application or computer program 1510 written in a programming language such as C, C++, Assembly, SQL, PYTHON, PROLOG, MATLAB, RUBY, RAILS, HASKELL, or other language to be translated into processor 1504 readable code. Alternatively, the compiler 1512 may be an interpreter that executes instructions/source code directly, translates source code into an intermediate representation that is executed, or that executes stored precompiled code. Such source code may be written in a variety of programming languages such as JAVA, JAVASCRIPT, PERL, BASIC, etc. After completion, the application or computer program 1510 accesses and manipulates data accepted from I/O devices and stored in the memory 1506 of the computer 1502 using the relationships and logic that were generated using the compiler 1512.

The computer 1502 also optionally comprises an external communication device such as a modem, satellite link, Ethernet card, or other device for accepting input from, and providing output to, other computers 1502.

In one embodiment, instructions implementing the operating system 1508, the computer program 1510, and the compiler 1512 are tangibly embodied in a non-transitory computer-readable medium, e.g., data storage device 1520, which could include one or more fixed or removable data storage devices, such as a zip drive, floppy disc drive 1524, hard drive, CD-ROM drive, tape drive, etc. Further, the operating system 1508 and the computer program 1510 are comprised of computer program 1510 instructions which, when accessed, read and executed by the computer 1502, cause the computer 1502 to perform the steps necessary to implement and/or use the present invention or to load the program of instructions into a memory 1506, thus creating a special purpose data structure causing the computer 1502 to operate as a specially programmed computer executing the protocol or method steps described herein. Computer program 1510 and/or operating instructions may also be tangibly embodied in memory 1506 and/or embodied in or coupled to source 1530 of the pulses 202 comprising electromagnetic fields (e.g., 1530 may comprise sources 500, 506), thereby making a computer program product or article of manufacture according to the invention. As such, the terms “article of manufacture,” “program storage device,” and “computer program product,” as used herein, are intended to encompass a computer program accessible from any computer readable device or media. Computer 1500 may comprise or be coupled to 1530.

Of course, those skilled in the art will recognize that any combination of the above components, or any number of different components, peripherals, and other devices, may be used with the computer 1502.

FIG. 16 schematically illustrates a typical distributed/cloud-based computer system 1600 using a network 1604 to connect client computers 1602 to server computers 1606. A typical combination of resources may include a network 1604 comprising the Internet, LANs (local area networks), WANs (wide area networks), SNA (systems network architecture) networks, or the like, clients 1602 that are personal computers or workstations (as set forth in FIG. 15 ), and servers 1606 that are personal computers, workstations, minicomputers, or mainframes (as set forth in FIG. 15 ). However, it may be noted that different networks such as a cellular network (e.g., GSM [global system for mobile communications] or otherwise), a satellite based network, or any other type of network may be used to connect clients 1602 and servers 1606 in accordance with embodiments of the invention.

A network 1604 such as the Internet connects clients 1602 to server computers 1606. Network 1604 may utilize ethernet, coaxial cable, wireless communications, radio frequency (RF), etc. to connect and provide the communication between clients 1602 and servers 1606. Further, in a cloud-based computing system, resources (e.g., storage, processors, applications, memory, infrastructure, etc.) in clients 1602 and server computers 1606 may be shared by clients 1602, server computers 1606, and users across one or more networks. Resources may be shared by multiple users and can be dynamically reallocated per demand. In this regard, cloud computing may be referred to as a model for enabling access to a shared pool of configurable computing resources.

Clients 1602 may execute a client application or web browser and communicate with server computers 1606 executing web servers 1610. Such a web browser is typically a program such as MICROSOFT INTERNET EXPLORER/EDGE, MOZILLA FIREFOX, OPERA, APPLE SAFARI, GOOGLE CHROME, etc. Further, the software executing on clients 1602 may be downloaded from server computer 1606 to client computers 1602 and installed as a plug-in or ACTIVEX control of a web browser. Accordingly, clients 1602 may utilize ACTIVEX components/component object model (COM) or distributed COM (DCOM) components to provide a user interface on a display of client 1602. The web server 1610 is typically a program such as MICROSOFT'S INTERNET INFORMATION SERVER.

Web server 1610 may host an Active Server Page (ASP) or Internet Server Application Programming Interface (ISAPI) application 1612, which may be executing scripts.

Generally, these components 1600-1616 all comprise logic and/or data that is embodied in/or retrievable from device, medium, signal, or carrier, e.g., a data storage device, a data communications device, a remote computer or device coupled to the computer via a network or via another data communications device, etc. Moreover, this logic and/or data, when read, executed, and/or interpreted, results in the steps necessary to implement and/or use the present invention being performed.

Although the terms “user computer”, “client computer”, and/or “server computer” are referred to herein, it is understood that such computers 1602 and 1606 may be interchangeable and may further include thin client devices with limited or full processing capabilities, portable devices such as cell phones, notebook computers, pocket computers, multi-touch devices, and/or any other devices with suitable processing, communication, and input/output capability.

Of course, those skilled in the art will recognize that any combination of the above components, or any number of different components, peripherals, and other devices, may be used with computers 1602 and 1606. Embodiments of the invention are implemented as a software protocol application on a client 1602 or server computer 1606. Further, as described above, the client 1602 or server computer 1606 may comprise a thin client device or a portable device that has a multi-touch-based display.

Process Steps

Method of Making a Register

FIG. 17A illustrates a method of making a system for implementing a quantum register. The method comprises the following steps.

Block 1700 represents obtaining or providing a device 1500 for coupling a qubit to a register. The device comprises one or more circuits or a computer 1502 configured to control a protocol 200 comprising a sequence 201 of pulses 202 synchronized with an RF field 204. Controlling the protocol comprises configuring (e.g., selecting, setting, or programming) a timing (e.g., spacing; σ/4 relative to other pukes and RF field), a phase (+/−x, +/−y), and a duration (pi, pi/2) of each of the pulses comprising a single qubit gate, a period τ, 210, and amplitude B^(RF) of the RF field, and a number of cycles M of the sequence, so that application of the protocol 200 controls a coherent spin exchange interaction {tilde over (Ŝ)}₊Î⁻+{tilde over (Ŝ)}⁻Î₊ between a register 206 and a qubit 208 having a zero magnetic dipole moment.

In one or more embodiments, the device comprises at least one of a signal generator, arbitrary waveform generator (e.g., comprising FPGA and digital to analog converter), or amplifier comprising the one or more circuits (e.g., as an embedded circuit or processor) outputting control signals that are used to control the output of the pulses (comprising the electromagnetic fields) and the RF field from one or more sources (e.g., lasers, microwave sources, or RF generator). In one or more examples, the sources of the pulses and RF field (e.g., the lasers) and microwave source(s) and RF source) comprise the one or more circuits, e.g., as an embedded system or processor, e.g., so as to form smart or programmable sources. The one or more circuits may be in central controller or distributed among the sources. In one or more examples, the arbitrary waveform generator (AWG) comprises the microwave sources and RF sources outputting the microwave pulses and RF field, and the AWG outputs the timing control signals to the laser sources.

In one or more examples, the device comprises a computer comprising or coupled to one or more processors; one or more memories; and one or more programs stored in the one or more memories, wherein the one or more programs executed by the one or more processors control the implementation of the protocol.

In one or more examples, the device comprises an application specific integrated circuit or field programmable gate array controlling the implementation of the protocol. In one or more examples, the one or more circuits comprise one or more timing circuits or a clock or a clock signal generator.

Block 1702 represents optionally coupling the device to one or more sources of electromagnetic fields. The one or more sources output the pulses comprising an electromagnetic field having a frequency (e.g., f_(g) in FIG. 5 a ) tuned to excite a transition between the first spin state and the second spin state.

Block 1704 represents optionally coupling the one or more sources to a photonic cavity.

Block 1706 represents optionally coupling the one or more sources to the qubit coupled to the register, e.g., via the photonic crystal. In one or more examples the qubit and the register are coupled, combined, or integrated with the photonic cavity.

The qubit comprises a first spin state (e.g., |0g>) having a zero magnetic dipole moment and a second spin state (e.g., |1g> having a zero magnetic dipole moment). The register comprises multiple register spins 100 having an energy level structure 102, wherein the register spins are indistinguishable so as to be configurable in basis states including a superposition state |Wv> used for storing the quantum state of the qubit.

A variety of systems including, but not limited to, solid state materials, can be used to implement the qubit and the register. In one example, the system comprises a spin carrying defect (e.g., an ion or nitrogen vacancy) in a host lattice (e.g., a crystal), wherein the spin carrying defect comprises the qubit and the host lattice comprises the register. Various rare earth doped crystals can be used. In one or more examples, the qubit ion comprising the qubit is Yb, Er, or Eu doped in a host crystal comprising register ions 122 surrounding the qubit ion. Examples include, but are not limited to, Yb:YVO (as described in the first example), Er:Y2SiO₅, or Eu:Y₂SiO₅). In another example, the system comprise a quantum dot in a host lattice, wherein the quantum dot (e.g., InGaAs or other semiconductor quantum dot) comprises the qubit and the host lattice (e.g., InGaAs or other semiconductor) comprises the register.

Block 1708 represents optionally coupling the qubit to a detector.

Block 1710 represents the end result, a system for coupling the qubit to a register. The system can be embodied in many ways, including, but not limited to, the following examples.

1. FIG. 15 , FIG. 2 , and FIG. 1 illustrate examples of a means for, or a device 1500 for coupling a qubit to a register, comprising a circuit or computer 1502 controlling a protocol 200 comprising a sequence 201 of pulses 202 synchronized with an RF field 204. Controlling the protocol comprises configuring (e.g., selecting, setting, or programming) a timing (e.g., spacing τ/4 relative to other pulses and RF field), a phase (+/−x., +/−y), and a duration (pi, pi/2) of each of the pulses comprising a single qubit gate, a period T and amplitude B^(RF) of the RF field, and a number of cycles M of the sequence, so that application of the protocol 200 controls a coherent spin exchange interaction (e.g., {tilde over (Ŝ)}₊Î⁻+{tilde over (Ŝ)}⁻Î₊) between a register 206 and a qubit 208 having a zero magnetic dipole moment: The qubit comprises a first spin state (e.g., |0g>) having a zero magnetic dipole moment and a second spin state (e.g., |1g> having a zero magnetic dipole moment. The register comprises multiple register spins 100 having an energy level structure 102, wherein the register spins are indistinguishable so as to be configurable in basis states including a superposition state |Wv> used for storing the quantum state of the qubit.

2. The device of example 1, wherein the protocol is configured to:

suppress or cancel one or more non-exchange interactions between the register and the qubit,

suppress or cancel noise coupled to the qubit and causing decoherence of a quantum state of the qubit,

enable the coherent spin exchange interaction that performs a quantum logic gate (e.g., a Clifford gate Ur as illustrated in FIG. 6 ), coherently transferring a quantum state of the qubit between the register and qubit.

The non-exchange interactions arise when the S_(x)I_(x) interaction is expressed in a form comprising spin preserving parts (spin exchange) and also non spin preserving parts (corresponding to the non-exchange interaction).

3. FIG. 17B illustrates an example wherein of the device of example 1 or 2, wherein the circuit controls:

application of a period of the protocol within a time period shorter than a rate of change of a magnetic noise (e.g., Overhauser field), so that the magnetic noise is quasistatic during the application of the period of protocol, the magnetic noise causing qubit decoherence and inducing a second order interaction (incoherent or random interaction) between the qubit and the register; and

at least one of a phase, duration, or time spacing of the pulses in the period so that:

-   -   one or more spin exchange interactions induced by the RE field         are preserved or maintained across the period;     -   one or more non-exchange interactions induced by the RE field         are cancelled across the period (e.g., components of the non         exchange interactions induced at different time instances in the         period cancel each other, or average to zero, over the period);     -   one or more (or any) exchange and one or more (or any)         non-exchange interactions induced by the magnetic noise are         cancelled across the period (e.g., components of these         interactions induced at different time instances in the period         cancel each other, or average to zero, over the period); and     -   the qubit decoherence induced by the magnetic noise is cancelled         over the period (e.g., decoherence induced at different time         instances in the period cancel each other, or average to zero,         over the period);     -   the RE field toggling between two values 214 of equal magnitude         and opposite polarity such that:         -   the period is associated with a frequency of a precession of             each of the multiple register spins about a predetermined             quantization axis (e.g., determined by an electric field             gradient generated by the host lattice at zero magnetic             field, or the application of a magnetic field); and         -   the amplitude is selected for a predetermined magnitude of             the coherent spin exchange interaction between the register             spins and qubit; and

so as to form a predictable (e.g., controllable, non random, deterministic) and coherent spin exchange interaction.

4. The device of any of the examples 1-3, wherein each of the single qubit gates comprises one of the pulses having the frequency (e.g., f_(g) in FIG. 5 a ) and duration (e.g., pi or pi/2) tuned to drive a transition between the first spin state and the second spin state.

6. The device of any of the examples claim 1-5 comprising a quantum memory 104 wherein the circuit:

controls application of the protocol in combination with are initialization of the qubit so as to configure the register spins in a polarized state |0v>;

controls application of one or more of the pulses to set a quantum state 106 of the qubit; and

controls application of the protocol so as to apply a first swap gate 108 (two qubit gate) transferring the quantum state of the qubit from the qubit to the register, thereby changing the polarized state to a corresponding state 110 of the register spins corresponding to the quantum state; and.

apply a second swap gate 112 retrieving the quantum state in the qubit from the register, thereby changing the corresponding state of the register spins to the polarized state.

6. The device of any of the examples 1-5, wherein configuring the register spins in the polarized state comprises polarizing the register, which is initially in an unpolarized state comprising any configuration of excitations of the register spins, by:

(a) initializing the qubit in the first spin state by controlling application of one or more initialization pulses of one or more initialization electromagnetic fields having one or more frequencies e.g., A, F, and fe in FIG. 5 a ) and tuned to initialize the quantum state of the qubit in the first spin state;

(d) applying the protocol transferring a spin excitation from the register spins to the qubit; and

(e) repeating steps (a) and (b) until all excitations of the register spins are transferred from the register to the qubit and the register spins are initialized in the polarized state, as characterized by a measurement of the qubit remaining in the first spin state after step (b).

The device of example 5 or 6, wherein the circuit controls application of the protocol so as to apply the first swap gate mapping (via the coherent spin exchange interaction) between the qubit and the register, such that:

if the qubit is in the first spin state, the corresponding state of the register is the polarized state |0v>,

if the qubit is in the second spin state, the corresponding state of the register is a W state |Wv>, and

if the qubit is in a superposition of the first spin state and the second spin state, the corresponding state of the register is a superposition 110 of the polarized state and the W state, and

wherein the W state is a superposition of all single spin excitation states of the register spins.

8. The device of any of the examples 1-7, wherein the circuit:

controls application of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state;

controls application of one or more of the pulses to set a quantum state of the qubit;

controls application of the protocol so as to apply a first square root of swap gate entangling the qubit with the register so as to form a Bell state; and

controls application of the protocol so as to apply a second square root of swap gate interacting with the Bell state so as to perform a measurement of the Bell state.

9. A repeater in a quantum network comprising the device of example 8.

10. FIG. 1 , FIG. 2 , and FIG. 5 illustrate an example of a system 112 comprising the device 1500 of any of the examples 1-9, further comprising:

a photonic cavity 114 coupled to a solid state material comprising the qubit and the register;

one or more microwave sources 500, 502 coupled to the qubit via a microwave waveguide, the microwave sources outputting one or more first microwave pulses fe and/or one or more second microwave pulses fg;

a radio frequency source 504 outputting the RE field; and

one or more laser sources 506, 508 outputting one or more laser pulses coupled to the qubit through the photonic cavity; and wherein:

the circuit controls the one or more laser sources and the one or more microwave sources so as to:

output initialization pulses comprising at least one of the one or more laser pulses A, F, or the one or more first microwave pulses fe having initialization frequencies for exciting one or more transitions initializing the qubit;

apply the protocol 200 comprising the single qubit gates comprising the second microwave pulses in synchronization with the RE field; and

output one or more readout electromagnetic fields having a readout frequency A for exciting a readout transition from the second spin state to a readout state |0e>, so as to stimulate output of readout pulses 116 from the readout state.

11. The device of any of the examples 1-10, wherein:

the pulses each comprise a pi pulse or a pi/2 pulse having at least one phase selected from +x, −x, +y, or −y, and

the circuit controls:

the sequence such that the period of the RF field is 2τ and a spacing of the pulses is τ/4, and

for a given magnitude of the spin exchange interaction determined by the amplitude of the RF field, a number of repeats M of the protocol that applies at least one of a swap gate transferring a quantum state between the qubit and the register, a square root of a swap gate for forming or measuring a Bell state, or that can be used to polarize the spins into a polarized state in combination with an initialization of the qubit.

12. The device of any of the examples 1-11 wherein the circuit selects and sets the duration and the timing of each of the pulses and a toggling of the RF field to engineer the coherent spin-exchange interaction comprising:

{tilde over (Ŝ)}₊Î⁻+{tilde over (Ŝ)}⁻Î₊,

where {tilde over (Î)}₊=|↑

↓|,{tilde over (Î)}⁻=|↓

↑| are the raising and lowering operators in an effective nuclear two-level manifold of the multiple spins in the register and {tilde over (Ŝ)}₊ are similarly defined for the qubit.

13. The device of any of the examples, wherein:

-   -   the RF field induces an interaction between the qubit and the         register comprising S_(z)I_(z) and at least one of S_(x)I_(x) or         S_(y)I_(y) including exchange and non-exchange components, where         Sx, Sy, Sx are the spin operators for the qubit ion and Ix, Iy,         Iz are the spin operators for the register ions along the x, y,         z cartesian axes respectively,     -   the control circuit applies the protocol that engineers the         interaction comprising only the coherent spin exchange         interaction by causing a cancelation of any non-exchange         components over the period, and     -   the pulses are synchronized with a precession of the register         ions about a predetermined quantization axis.

14. FIG. 2 illustrates an example of the device of any of the examples 1-13, wherein the RF field comprises a square wave and the sequence of pulses comprise:

in a first half period τ of the square wave a sequence of the second pulses comprising:

a first pi/2 pulse having a phase +Y followed by a first pi pulse having a phase +Y, the beginning of the first pi/2 pulse and the center of the first pi pulse separated in time by τ/4;

a second pi/2 pulse immediately followed by a third pi/2 pulse, the end of the second pi/2 pulse separated in time from the center of the first pi pulse by τ/4, wherein the second pi/2 pulse has a phase −Y and the third pi/2 pulse has a phase −X;

a second pi pulse having a phase −X and following the third pi/2 pulse, a center of the second pi puke separated in time from the center of the first pi pulse by τ/2; and

a fourth pi/2 pulse having a phase −X, wherein the end of the fourth pi/2 pulse is separated in time from center of the second pi pulse by τ/4; and

in a second half period τ of the square wave, a repeat of the sequence of second pulses but wherein the first pi/2 pulse, the first pi pulse, and the second pi/2 pulse have opposite phase as compared to the first pi/2 pulse, the first pi pulse, and the second pi/2 puke in the first half period, respectively.

15. The device of any of the examples 1-14, wherein the protocol de-couples the qubit from decoherence noise and random interactions caused by a nuclear Overhauser field generated by a host lattice 118 in which the qubit is located.

16. A system 112 for implementing a quantum register comprising the device of any of the examples claim 1-15 coupled to:

a spin carrying defect 120 in a host lattice, wherein the spin carrying defect comprises the qubit and the host lattice 118 comprises the register, or

a quantum dot in a host lattice, wherein the quantum dot comprises the qubit and the host lattice comprises the register.

17. The system of claim 16, wherein the spin carrying defect is a qubit ion comprising the qubit and the register comprises a lattice 118 of register ions 122 surrounding the qubit

18. The device of any of the examples 1-17, wherein the multiple register spins 100 in the register comprise nuclear spins and the first spin state and the second spin state comprise electron spin states.

19. The device of any of the examples, wherein the protocol controls oscillations between a first system state |1g>|0v>, representing a spin excitation in the qubit and register ions in the polarized state, and a second system state |0g>|W> where |W> is an entangled |W> spin state of the register comprising the spin excitation transferred from the qubit.

20. A protocol 200 for controlling a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment, wherein the qubit comprises a first spin state and a second spin state both having zero magnetic dipole moment; and the register comprises multiple indistinguishable spins. The protocol comprises a toggling RF field or magnetic field synchronized to a sequence of pulses, wherein a period (e.g. 2T) of the toggling RF field or magnetic field is matched to a spacing (e.g., T/4) of the pulses comprising single qubit gates (e.g., clifford gates performing unitary operations) and the protocol modulates the spin exchange interaction so as to transfer quantum information to or from the qubit.

21. The protocol of any of the examples, wherein the spin exchange interaction comprises an interaction between an electron spin of the qubit and a nuclear spin of the register (e.g., electron-nuclear dipole interaction) or an interaction between electron spins of the qubit and the register. Block 1712 represents optionally coupling the system in or to an application, e.g., in or to a quantum computer, in or to a quantum network, or in repeater for a quantum network.

22. The protocol of any of the examples 1-21, wherein the RE field comprises or is a magnetic field or the radio frequency (RF) field has a frequency in range 20 kHz-300 GHz.

Method of Performing: Qubit Operations with a Controlled Spin Exchange Interaction

FIG. 18 is a flowchart illustrating a method for coupling a qubit to a quantum register.

Block 1800 represents obtaining a protocol comprising a sequence of pukes synchronized with an RF field, the protocol further comprising a timing, a phase, and a duration of each of the pukes comprising a single qubit gate, and a period and amplitude of the RF field, wherein application of the protocol controls a coherent spin exchange interaction between a register and a qubit.

Block 1802 represents applying one or more cycles of the protocol to the qubit, so as to modulate the coherent spin exchange interaction transferring a spin excitation between the qubit and the register. The qubit comprises a first spin state and a second spin state both having a zero magnetic dipole moment, the register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing a quantum state of the qubit; and the pulses comprise an electromagnetic field tuned to excite a transition between the first spin state and the second spin state.

1. Quantum Memory

FIG. 19 illustrates a method of applying a number of cycles of the protocol so as transfer quantum information between the qubit and the register.

Block 1900 represents applying a first number of the cycles of the protocol to the qubit in combination with an initialization of the qubit so as to configure the register spins in a polarized state.

Block 1902 represents applying one or more of the pulses to the qubit to set a quantum state of the qubit.

Block 1904 represents applying a second number of the cycles of the protocol to the qubit so as to apply a first swap gate (two qubit gate) transferring a quantum state of the qubit from the qubit to the register, thereby changing the polarized state to a corresponding state of the register spins corresponding to the quantum state.

Block 1906 represents applying one or more cycles of the protocol to the qubit so as to apply a second swap gate retrieving the quantum state in the qubit from the register, thereby changing the corresponding state of the register spins to the polarized state.

2. Bell State measurement

FIG. 20 is a flowchart illustrating a method of forming and measuring Bell states. The method comprises the following steps.

Block 2000 represents applying a first number of the cycles of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized w state.

Block 2002 represents applying one or more of the pulses to the qubit to set a quantum state of the qubit.

Block 2004 represents applying a second number of the cycles of the protocol to the qubit so as to apply a first square root of swap gate entangling the qubit with the register so as to form a Bell state.

Block 2006 represents applying one or more cycles of the protocol to the qubit so as to apply a second square root of swap gate interacting with the Bell state so as to perform a measurement of the Bell state.

REFERENCES FOR SUPPLEMENTARY EXAMPLES

The following references are incorporated by reference herein.

-   [1] Kindem, J. M. et al. Characterization of Yb 3+171:YVO4 for     photonic quantum technologies. Phys. Rev. B 98, 1-10 (2018) -   [2] Bleaney, B., Gregg, J. F., De Oliveira, A. C. & Wells, M. R.     Nuclear magnetic resonance of 51 Vt(I=7/2) in lanthanide     vanadates: II. The nuclear electric quadrupole interaction. J.     phys., C, Solid state phys. 15, 5293-5303 (1982). -   [3] Bleaney, B., Gregg, J. F., De Oliveira, A. C. & Wells, M. R.     Nuclear magnetic resonance of 51 V (1=7/2) in lanthanide     vanadates: I. The paramagnetic shifts. J. phys., C, Solid state     phys. 15, 5293-5303 (1982). -   [4] Cohen-Tannoudji, C., Dupont-Roc, J. & Grynberg, G. Atom-Photon     Interactions (Wiley-VCH, Weinheim, 2004). -   [5] Bermudez, A., Jelezko, F., Plenio, M. B. & Retzker, A.     Electron-mediated nuclear-spin interactions between distant     nitrogen-vacancy centers. Phys. Rev. Lett. 107, 3-7 (2011). -   [6] Knill, E. et al. Randomized benchmarking of quantum gates.     Physical Review A—Atomic, Molecular, and Optical Physics 77, 1-7     (2008) -   [7] Gullion, T., Baker, D. B. & Conradi, M. S. New, compensated     Can-Purcell sequences. J. Magn. Reson. 89, 479-484 (1990) -   [8] Kindem, J. M. et al. Control and single-shot readout of an ion     embedded in a nanophotonic cavity. Nature 580, 201-204 (2020) -   [9] Taylor, J. M., Marcus, C. M. & Lukin, M. D. Long-Lived Memory     for Mesoscopic Quantum Bits. Phys. Rev. Lett. 90, 4 (2003) -   [10] Hartmann, S. R. & Hahn, E. L. Nuclear double resonance in the     rotating frame. Phys. Rev. 128, 2042-2053 (1962). -   [11] Scheuer, J. et al. Robust techniques for polarization and     detection of nuclear spin ensembles. Phys. Rev. B 96, 1-10 (2017) -   [12] Urbaszek, B. et al. Nuclear spin physics in quantum dots: An     optical investigation. Rev. Mod. Phys. 85, 79-133 (2013). -   [13] Bernien, H. et al. Heralded entanglement between solid-state     qubits separated by three metres. Nature 497, 86-90 (2013). -   [14] Further information on one or more embodiments of the present     invention can be found in     https://www.nature.com/articles/s41586-021-04293-6, Ruskuc, A., Wu,     C J., Rochman, J. et al. Nuclear spin-wave quantum register for a     solid-state qubit. Nature 602, 408-413 (2022).     https://doi.org/10.1038/s41586-021-04293-6. -   [15] US Patent Application Publication No. 20210028863 by Faraon et,     al., entitled Optical Quantum Networks with Rare Earth Ions. U.S.     patent application Ser. No. 16/937,379.

Advantages and Improvements

Usually, working with dense nuclear spin ensembles leads to qubits with short coherence times where information cannot be transferred or stored for long. Additionally, these nuclear spins are often indistinguishable, meaning information cannot be stored on a single nuclear spin (as is commonly a requirement with other systems/protocols).

Example systems described herein resolve the issue of maintaining qubit coherence by using a transition with no magnetic dipole moment. However, the lack of magnetic dipole moment also inhibits the interactions needed to transfer quantum information to the nuclear spins. The pulse sequence disclosed herein enables this interaction despite the lack of magnetic dipole moment. Additionally, the form of interaction (spin preserving) enables storage of information in a delocalized fashion across multiple indistinguishable nuclear spins (so that single nuclear spin storage is no longer a requirement).

Thus, advantages of the protocol disclosed herein include:

-   -   Enabling initialization and control of a multi-level nuclear         spin ensemble, which provides a much larger Hilbert space for         quantum simulation compared to conventional single spin-1/2         nuclei.     -   Providing a novel configuration of pulse sequences enabling         coherent control of the nuclear spin register using magnetically         insensitive (and hence low-noise) qubit transitions.     -   Enabling the realisation of a reproducible and deterministic         quantum register, a critical requirement for building scalable         quantum networks.

Definitions

As known to a person skilled in the art, a pi pulse may refer to a pulse of light (e.g., laser) or microwaves generally resonant with a transition between two levels, the pulse being calibrated via known methods to move the population/excitation fully from one level to another.

Accordingly, an optical pi pulse is a .pi. pulse in the optical (e.g., visible) domain/frequencies, and a microwave pi pulse is a .pi. pulse in the microwave domain/frequencies. It should be noted that a pi pulse can move (transfer) population/excitations with a probability of 1, so as to change the state of the qubit between the two spin states 0g and 1g, whereas as a non-pi pulse can transfer population/excitations with some probability between 0 and 1, and not necessarily 1, so as to form the qubit comprising a superposition of the spin states 0g and 1g.

In one or more examples, a spin-exchange interaction preserves total angular momentum of the system but may allow other aspects of the system to change. When two spins in the qubit and register experience a spin-exchange interaction, the total spin of the qubit-register system is preserved yet the orientation of the individual spins in the register and qubit may change. For example, if qubit A and register B are in opposite spin states, a spin-exchange interaction reverses the spins

A(↑)+B(↓)→A(↓)+B(↑)

CONCLUSION

This concludes the description of the preferred embodiment of the present invention. The foregoing description of one or more embodiments of the invention has been presented for the purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed. Many modifications and variations are possible in light of the above teaching. It is intended that the scope of the invention be limited not by this detailed description, but rather by the claims appended hereto. 

What is claimed is:
 1. A device for coupling a qubit to a register, comprising: a circuit for controlling application of one or more cycles of a protocol comprising a sequence of pulses synchronized with an FT field, a timing, a phase, and a duration of each of the pulses, and a period and amplitude of the RF field, wherein: application of the one or more cycles of the protocol controls a coherent spin exchange interaction between a register and a qubit having a zero magnetic dipole moment; the qubit comprises a first spin state and a second spin state both of which have the zero magnetic dipole moment; the register comprises multiple register spins having an energy level structure, the register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing the quantum state of the qubit, and the pulses each comprise an electromagnetic field tuned to excite a transition between the first spin state and the second spin state.
 2. The device of claim 1, wherein the protocol is configured to: suppress or cancel one or more non-exchange interactions between the register and the qubit, suppress or cancel noise coupled to the qubit and causing decoherence of a quantum state of the qubit, enable the coherent spin exchange interaction that performs a quantum logic gate, coherently transferring a quantum state of the qubit between the register and qubit.
 3. The device of claim 1, wherein the circuit controls: application of a period of the protocol within a time period shorter than a rate of change of a magnetic noise, so that the magnetic noise is quasistatic during the application of the period of the protocol, wherein the magnetic noise causes qubit decoherence and induces a second order interaction (incoherent interaction) between the qubit and the register; and at least one of a phase, duration, or time spacing of the pulses in the period so that: the spin-exchange interactions induced by the RF field are preserved or maintained across the period; one or more non-exchange interactions induced by the RF field are cancelled across the period; any exchange interactions and any non-exchange interactions induced by the magnetic noise are cancelled across the period; and the qubit decoherence induced by the magnetic noise is cancelled over the period; and the RF field toggling between two values of equal magnitude and opposite polarity such that: the period is associated with a frequency of a precession of each of the multiple register spins about a predetermined quantization axis; and the amplitude is selected for a predetermined magnitude of the coherent spin exchange interaction between the register spins and qubit; and so as to form the coherent spin exchange interaction.
 4. The device of claim 3, wherein the protocol comprises the sequence of single qubit gates, each of the single qubit gates comprising one of the pulses having the frequency and duration tuned to drive a transition between the first spin state and the second spin state.
 5. The device of claim 1 comprising a quantum memory, wherein the circuit: controls application of a number of cycles of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state; controls application of one or more of the pulses to set a quantum state of the qubit; and controls application of a number of cycles of the protocol so as to apply a first swap gate (two qubit gate) transferring the quantum state of the qubit from the qubit to the register, thereby changing the polarized state to a corresponding state of the register spins corresponding to the quantum state; and controls application of a number of cycles of the protocol so as to apply a second swap gate retrieving the quantum state in the qubit from the register, thereby changing the corresponding state of the register spins to the polarized state.
 6. The device of claim 5, wherein configuring the register spins in the polarized state comprises polarizing the register, which is initially in an unpolarized state comprising any configuration of excitations of the register spins, by: (a) initializing the qubit in the first spin state by controlling application of one or more initialization pulses of an initialization electromagnetic field having a frequency tuned to initialize the quantum state of the qubit in the first spin state; (b) applying one or more cycles of the protocol transferring a spin excitation from the register spins to the qubit; and (c) repeating steps (a) and (b) until all excitations of the register spins are transferred from the register to the qubit and the register spins are initialized in the polarized state, as characterized by a measurement of the qubit remaining in the first spin state after step (b).
 7. The device of claim 5, wherein the circuit controls application of the protocol so as to apply the first swap gate mapping (via the coherent spin exchange interaction) between the qubit and the register, such that: if the qubit is in the first spin state, the corresponding state of the register is the polarized state, if the qubit is in the second spin state, the corresponding state of the register is a W state, and if the qubit is in a superposition of the first spin state and the second spin state, the corresponding state of the register is a superposition of the polarized state and the W state, and wherein the W state is a superposition of all single spin excitation states of the register spins.
 8. The device of claim 1, wherein the circuit: controls application of one or more cycles of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state; controls application of one or more of the pulses to set a quantum state of the qubit; controls application of one or more cycles of the protocol so as to apply a first square root of swap gate entangling the qubit with the register so as to form a Bell state; and controls application of one or more cycles of the protocol so as to apply a second square root of swap gate interacting with the Bell state so as to perform a measurement of the Bell state.
 9. A repeater in a quantum network comprising the device of claim
 8. 10. The repeater of claim 9, further comprising: a photonic cavity coupled to a solid state material comprising the qubit and the register; one or more microwave sources coupled to the qubit via a microwave waveguide, the microwave sources outputting one or more first microwave pulses and/or one or more second microwave pulses; a radio frequency source outputting the RF field; and one or more laser sources outputting one or more laser pulses coupled to the qubit through the photonic cavity; and wherein: the circuit controls the one or more laser sources and the one or more microwave sources so as to: output initialization pulses comprising at least one of the one or more laser pulses or the one or more first microwave pulses having initialization frequencies for exciting one or more transitions initializing the qubit; apply the protocol comprising the single qubit gates comprising the second microwave pulses in synchronization with the RF field; and output one or more readout electromagnetic fields having a readout frequency for exciting a readout transition from the second spin state to a readout state, so as to stimulate output of third pulses from the readout state.
 11. The device of claim 1, wherein: the pulses each comprise a pi pulse or a pi/2 pulse having at least one phase selected from +x, −x, +y, or −y, and the circuit controls: the sequence such that the period of the RF field is 2τ, and a spacing of the pulses is τ/4, and for a given magnitude of the spin exchange interaction determined by the amplitude of the RF field, a number of repeats or cycles of the protocol that applies at least one of a swap gate transferring a quantum state between the qubit and the register, a square root of a swap gate for forming or measuring a Bell state, or that can be used to polarize the spins into a polarized state in combination with an initialization of the qubit.
 12. The device of claim 1 wherein the circuit selects the duration and the timing of each of the pulses and a toggling of the RF field to engineer the coherent spin-exchange interaction comprising: {tilde over (Ŝ)}₊Î⁻+{tilde over (Ŝ)}⁻Î₊, where {tilde over ({right arrow over (I)})}+=|↑

↓|,{acute over (Î)}⁻=|↓

↑| are the raising and lowering operators in an effective nuclear two-level manifold of the multiple spins in the register and {tilde over (Ŝ)}₊ are similarly defined for the qubit.
 13. The device of claim 1, wherein the RF field comprises a square wave and the sequence of pulses comprise: in a first half period T of the square wave a sequence of the second pulses comprising: a first pi/2 pulse having a phase +Y followed by a first pi pulse having a phase +Y, the beginning of the first pi/2 pulse and the center of the first pi pulse separated in time by τ/4; a second pi/2 pulse immediately followed by a third pi/2 pulse, the end of the second pi/2 pulse separated in time from the center of the first pi pulse by τ/4, wherein the second pi/2 pulse has a phase −Y and the third pi/2 pulse has a phase −X; a second pi pulse having a phase −X and following the third pi/2 pulse, a center of the second pi pulse separated in time from the center of the first pi pulse by τ/2; and a fourth pi/2 pulse having a phase −X, wherein the end of the fourth pi/2 pulse is separated in time from center of the second pi pulse by τ/4; and in a second half period τ of the square wave, a repeat of the sequence of second pulses but wherein the first pi/2 pulse, the first pi pulse, and the second pi/2 pulse have opposite phase as compared to the first pi/2 pulse, the first pi pulse, and the second pi/2 pulse in the first half period, respectively.
 14. The device of claim 1, wherein the protocol de-couples the qubit from decoherence noise and random interactions caused by a nuclear Overhauser field generated by a host lattice in which the qubit is located.
 15. A system for implementing a quantum register comprising the device of claim 1 coupled to: a spin carrying defect in a host lattice, wherein the spin carrying defect comprises the qubit and the host lattice comprises the register, or a quantum dot in a host lattice, wherein the quantum dot comprises the qubit and the host lattice comprises the register.
 16. The system of claim 15, wherein the spin carving defect is a qubit ion comprising the qubit and the register comprises a lattice of register ions surrounding the qubit ion.
 17. The device of claim 1, wherein the multiple spins in the register comprise nuclear spins and the first spin state and the second spin state comprise electron spin states.
 18. A method for coupling a qubit to a quantum register, comprising: obtaining a protocol comprising a sequence of pulses synchronized with an RF field, the protocol further comprising a timing, a phase, and a duration of each of the pulses comprising a single qubit gate, and a period and amplitude of the RF field, wherein application of the protocol controls a coherent spin exchange interaction between a register and a qubit; and applying one or more cycles of the protocol to the qubit, so as to modulate the coherent spin exchange interaction transferring a spin excitation between the qubit and the register; and wherein: the qubit comprises a first spin state and a second spin state both having a zero magnetic dipole moment, the register spins are indistinguishable so as to be configurable in basis states including a superposition state used for storing a quantum state of the qubit; and the pulses comprise an electromagnetic field tuned to excite a transition between e first spin state and the second spin state.
 19. The method of claim 18, further comprising applying a number of cycles of the protocol so as transfer quantum information between the qubit and the register, comprising: applying a first number of the cycles of the protocol to the qubit in combination with an initialization of the qubit so as to configure the register spins in a polarized state; applying one or more of the pulses to the qubit to set a quantum state of the qubit; applying a second number of the cycles of the protocol to the qubit so as to apply a first swap gate (two qubit gate) transferring a quantum state of the qubit from the qubit to the register, thereby changing the polarized state to a corresponding state of the register spins corresponding to the quantum state; and applying one or more cycles of the protocol to the qubit so as to apply a second swap gate retrieving the quantum state in the qubit from the register, thereby changing the corresponding state of the register spins to the polarized state.
 20. The method of claim 18, further comprising a number of the cycles of the protocol so as to form and measurement of a Bell state, comprising: applying a first number of the cycles of the protocol in combination with an initialization of the qubit so as to configure the register spins in a polarized state; applying one or more of the pulses to the qubit to set a quantum state of the qubit; applying a second number of the cycles of the protocol to the qubit so as to apply a first square root of swap gate entangling the qubit with the register so as to form a Bell state; and applying one or more cycles of the protocol to the qubit so as to apply a second square root of swap gate interacting with the Bell state so as to perform a measurement of the Bell state. 